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ABSTRACT 

We assemble a sample of 24 hydrogen-poor super-luminous supernovae (SLSNe). Pa¬ 
rameterizing the light curve shape through rise and decline timescales shows that 
the two are highly correlated. Magnetar-powered models can reproduce the correla¬ 
tion, with the diversity in rise and decline rates driven by the diffusion timescale. 
Circumstellar interaction models can exhibit a similar rise-decline relation, but only 
for a narrow range of densities, which may be problematic for these models. We find 
that SLSNe are approximately 3.5 magnitudes brighter and have light curves 3 times 
broader than SNe Ibc, but that the intrinsic shapes are similar. There are a num¬ 
ber of SLSNe with particularly broad light curves, possibly indicating two progenitor 
channels, but statistical tests do not cleanly separate two populations. The general 
spectral evolution is also presented. Velocities measured from Fe II are similar for 
SLSNe and SNe Ibc, suggesting that diffusion time differences are dominated by mass 
or opacity. Flat velocity evolution in most SLSNe suggests a dense shell of ejecta. If 
opacities in SLSNe are similar to other SNe Ibc, the average ejected mass is higher by 
a factor 2-3. Assuming k = 0.1 cm 2 g -1 , we estimate a mean (median) SLSN ejecta 
mass of 10 M 0 (6 M 0 ), with a range of 3-30 M 0 . Doubling the assumed opacity brings 
the masses closer to normal SNe Ibc, but with a high-mass tail. The most proba¬ 
ble mechanism for generating SLSNe seems to be the core-collapse of a very massive 
hydrogen-poor star, forming a millisecond magnetar. 

Key words: Supernovae: general - Supernovae: LSQ14mo - Supernovae: LSQ14bdq 
- Supernovae: SN 2013hx 


1 INTRODUCTION 

By the end of the last century, the diversity in observed 
supernovae (SNe) could be explained in terms of well- 
understood physical differences: explosions due to iron core 
collapse, or thermonuclear runaway; progenitor stars being 
hydrogen-rich, or -poor; the presence or absence of a dense 


circumstellar medium (e.g. Filippenko 1997). However in 


the era of expansive, untargetted sky surveys, such as the 


Palomar Transient Factory (PTF; 

Rau et al. 2009) Pan- 

STARRSl (PS1; Kaiser et al. 2010 

), the La Silla QUEST 


Time Transient Survey (CRTS; |Drake et al.|2009 ) thousands 


of supernovae are being discovered each year, revealing un¬ 
expected new types of explosions. 

Perhaps the most mysterious of these are the “super- 
luminous” supernovae (SLSNe; |Gal- Yam 2012), so called be¬ 


cause of peak luminosities over 2 magnitudes brighter than 
the bulk of the SN population (of which they make up only 
ffi.01%; e.g. Quimby et al.||2013 McCrum et al.||2015 ). Al¬ 


though outliers had been noted previously, | Quimby et a h] 


(2011) were the first to define the SLSN class, and showed 


that the light curves of these objects are difficult to explain 
using only the release of energy deposited by the SN shock 
wave and the decay of synthesised 56 Ni, which are the energy 
sources for the known SN population. Because of their high 
lumosities, SLSNe have been spectroscopically confirmed up 
to redshift z > 1.5 (Be rger et al.|2012 ), and candidates have 


even been detected photometrically at z — 2 — 4 (Cooke 
et al.|20l2 ) 


Three main models have been proposed to account for 
the enormous energy radiated by SLSNe. One is a central 
engine, such as a magnetised neutron star spinning with a 
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period of order milliseconds (often referred to as the magne¬ 
tar model; Kasen &; Bildsten]|2010| |Woosley||2010 ). Magne- 
tars (albeit with longer periods) have been observed in our 
Galaxy, and are thought to originate from stars with main- 
sequence masses Mzams = 30 — 40 M 0 ( Gaensler et al.|2 005 


but see also Davies et al. 2009). The compact object spins 


down and heats the ejected SN gas through high-energy 
emission, though it is an open question as to exactly how this 
energy can thermalise in the ejecta (Metzger & Piro 2014). 


Another model is a collision with (or shock breakout from) a 
highly opaque circumstellar medium (CSM), releasing shock 
energy at a large radius (jWoosle y, Blinni kov &; Heger||2007| 


Ofek et al.||2010| | Chevalier &; Irwin||2011| |Ginzburg &; Bal-| 

berg|2012] k Comparisons with data have shown that magne¬ 

tar engines and circumstellar interaction can both reproduce 
the range of shapes of SLSN light curves and distinguishing 


between them has been problematic ([inserra et al. 2013 
Chatzopoulos et al.||2013 Nicholl et al.|2014). 


The final mechanism is thermonuclear runaway in a star 
with a main-sequence mass above 13OM 0 , triggered by pair- 
production in the hot carbon-oxygen core (a pair-instability 
supernova, or PISN; Barkat, Rakavy & Sack 1967; Rakavy & 
Shaviv|1967 ). A small number of SLSNe have been proposed 


to be PISNe based on their slowly fading light curves, which 
have decline rates tha t approximately match the decay of 


s Co (e.g. SN 2007bi; Gal-Yam et al. 


2009 


Young et al. 


2010| ). This physical scenario may be responsible for only 


a fraction of the SLSN population, since the slowly fading 
types appear to be rarer than those which decay too rapidly 


to be radioactively powered (Nicholl et al. 2013 McCrum 


|et al. ||20151 estimate the slowly fading SLSNe to be around 
10% of the total SLSN population). However the physical re¬ 
ality of pair-instability explosions is not firmly established, 
and the rise times of well-observed events do not satisfac¬ 
torily match the predictions of quantitative light curve and 
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spectral modelling (Nichollet al.[ 2013; Mc Crum et al.|2014 up. The first studies of SLSN samples were recently con- 


Dessart et al.||2012 ). 

Extensive studies of SLSNe in the last few years il¬ 
lustrate a diversity of spectral and photometric proper¬ 


ties. Most display no signs of hydrogen ( Pastorello et al. 
Inserra et al.|2013 ) and occur in low-metallicity dwarf 


2010 


galaxies - perhaps similar to the hosts of long-duration 
gamma-ray bursts ( Neill et al.|2011 Chen et al.|2013 Lun- 


nan et al. 2014), and possibly in even more extreme star¬ 


forming environments (Chen et al. 2014 Leloudas et al. 


2015). These objects have been termed SLSNe Ic, by anal¬ 
ogy with normal-luminosity SNe from stripped progenitors. 
Within this group, objects can show nearly identical spectral 
evolution and yet have very different light curves (Inserra 
et al.|2013| [Nicholl et al.|2014l ). Some SLSNe Ic (S N 2007bi- 
like objects) clearly evolve on very long timescales (| Gal-Yam 


et al. 1 2 009| | Nicholl et al.||2013 McCrum et al.|2014| ), but it 

is unclear whether these form a distinct subclass (and arise 
from a different physical mechanism, such as the pair insta¬ 
bility) or are part of a continuous distribution. For example, 
Inserra V Smartt (2014) found that their standardisation of 


SLSN peak magnitudes could encompass the slowly fading 
objects as well as the more typical ones. 

Another group, SLSNe II, do have hydrogen in their 
spectra. This is sometimes in the form of strong, multi- 
component emission lines, almost certainly indicating in¬ 
teraction with CSM. The prototypical example here is SN 
2006gy ( Ofek et al.|2007 Smith et al.|2007| ), and the physi¬ 
cal mechanism (the conversion of kinetic energy to radiative 
energy by shocks at the ejecta-CSM collision) is well es¬ 
tablished. Quite a few such objects are now known - e.g. 
SNe 2006tf ( Smith et al.||2008 ), 2008fz ( Drake et al.||2010 ), 
2008am ( Chatzopoulos et al.|2011 ) and 2003ma (Rest et al. 
2011) - and these are quite correctly dubbed ‘SLSNe Iln’, 


by analogy with the fainter SNe Iln, which are hydrogen- 
rich SNe showing narrow spectral lines from shocked CSM. 
However, a few SLSNe have much weaker hydrogen lines 
visible, which are not obviously multi-component and do 
not unambiguously point to interaction being the dominant 
power source of the radiative energy. The earliest example 


of this class is SN 2008es (Gezari et al. 2009 Miller et al. 


2009). Although they are classed as SLSNe II, they resem¬ 


ble SLSNe Ic with H lines superimposed. Their overall light 
curve and other spectral properties are closer to SLSNe Ic 
than to the SLSNe Iln. One object in particular, CSS121015, 


prompted Benetti et al. (2014) to propose that the two spec¬ 


troscopic classes of SLSNe may in fact come from the same 
underlying physical process, with their observational prop¬ 
erties modified by the hydrogen mass in the ejecta/CSM. 
CSS121015 had one of the highest peak luminosities of any 
SLSN to date, but its light curve could be well fit with both 
interaction and magnetar models. The authors favoured the 
interaction scenario, due to the presence of time-variability 
in the narrow Balmer emission, indicating slow-moving ma¬ 
terial close to the SN. 

As SLSNe are extremely rare events, and were largely 


unknown before the discoveries of SNe 2005ap (Quimby 
et al.|2007 ) and SCP06F6 ( Barbary et al.|2009 ), the paucity 
of observed events has so far restricted the analysis of their 
properties as a group. This situation is now beginning to 
change, as transient surveys are becoming better at picking 
out these objects and dedicating resources to follow them 


ducted by Inserra et al. (2013); Inserra V Smartt (2014), who 


showed that their properties could be explained by magnetar 
powered models, and their potential utility as standardisable 
candles for cosmology. 

The Public ESO Spectroscopic Survey of Transient Ob¬ 
jects (PESSTO) has a strategy designed to classify un¬ 
usual types of transients early in their evolution, using 
light curve information from feeder surveys such as La Silla 
QUEST (LSQ; |Baltay et al. 2013), Catalina Real-Time 
Transient Survey (CRTS; |Drake et al. 2009), SkyMapper 
(Keller et al. |2007| , OGLE-IV ( Wyrzykowski et al.[ 2014) 
and Pan-STARRSl (e.g. as described in Inserra et al.|2013 ). 
PESSTO is described in detail in Smartt et al. ~( 2014) and re¬ 


duced and calibrated spectra are publicly available through 
both ESO archiv^] and WISeREFp] (Weizmann Interactive 
Supernova data REPository; Yaron V Gal-Yam 2012). In 


this work, we construct the largest SLSN sample to date: a 


total of 24 objects, composed of 7 from PESSTO (Nicholl 


|et al.|2014||Benetti et al.|2014[ Inserra et al. in prep., Nicholl 
et al. in prep., Chen et al. in prep.), and 17 from the liter¬ 
ature. As all of the theoretical scenarios invoked to power 
SLSNe likely require stars more massive than typical SN 
progenitors, we investigate whether our sample have sys¬ 
tematically different ejected mass than normal-luminosity 
stripped-envelope SNe. In section [2] we describe the SLSNe 
in our sample, including the new PESSTO objects. The con¬ 
struction of our bolometric light curves is outlined in section 
[3] We investigate the light curve timescales in section]!] lead¬ 
ing to an analysis of generalised light curve shapes in section 
[5j and a search for evidence of a bimodal SLSN Ic popula¬ 
tion (i.e. with rapid and slow decline rates after maximum 
luminosity) in section [6| The typical spectral evolution is in¬ 
vestigated in section [7] Velocity measurements from spectra 
are described in section [8] and these are used to estimate 
SLSN masses relative to normal hydrogen-poor SNe Ibc in 
section [9] We summarise our main results in section [To] and 
conclude in section m 


2 THE SAMPLE 

In this work we focus on the Type Ic SLSNe. However, we 
also include three SLSNe II. While SLSNe Type Iln, such as 
SN 2006gy, show prominent multicomponent Balmer lines 
indicating circumstellar interaction, the three objects used 
in our sample showed only weak and/or broad hydrogen 
lines. As the power source for these objects is ambiguous, 
they may be related to SLSNe Ic. A full summary of the 
sample is given in Table ^ A comparison sample of normal- 
luminosity stripped-envelope SNe (Types lb, Ic and broad- 
lined Ic) is listed in Table [ 2 ] This contains a compilation 
of SNe Ibc from the literature that have good photometric 
data in griz , as well as the homogeneous SDSS II sample of 


Taddia et al. (2015) 


1 For details on how to get the PESSTO Phase 3 data, see 
www.pessto.org 

2 http://www.weizmann.ac.il/astrophysics/wiserep/ 
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Table 1 . SLSNe in our sample 


Name 

Type 

z 

M ■ * 

IVlgriZ 

Reference 



‘Gold’ sample: rest-frame gri(z) coverage 

SN2007bi 

let 

0.127 

-20.20 

Gal-Yam et al. (2009) 

SN2008es 

II 

0.205 

-21.43 

Gezari et al. ( 

2009) 

2009) 






Miller et al. L 


SN2010gx 

Ic 

0.230 

-20.64 

Pastorello et al.l 





(2010 

), 







Quim 

yy et al. 

(2011) 

SN2011ke 

Ic 

0.143 

-20.69 

Inserra et al. i 

2013 


SN2011kf 

Ic 

0.245 

-20.80 

Inserra et al. i 

2013 


SN2012il 

Ic 

0.175 

-20.73 

Inserra et al. i 

2013 


SN2013dg 

Ic 

0.265 

-20.30 

Nicholl et al. i 

20ll 

1 

SN2013hx 

II 

0.130 

-20.84 

Inserra et al. i 

(in 






prep) 




LSQ12dlf 

Ic 

0.255 

-20.68 

Nicholl et al. i 

2014) 

LSQ14mo 

Ic 

0.253 

-19.95 

Chen et al. (in prep) 

LSQ14bdq 

Ic 

0.347 

-21.68 

Nicholl et al. i 

2015 

1 

PTFlOhgi 

Ic 

0.100 

-19.61 

Inserra et al. i 

2013 


PTFllrks 

Ic 

0.190 

-20.01 

Inserra et al. i 

2013 


PTF12dam 

let 

0.107 

-20.56 

Nicholl et al. i 

2013 

1 

CSS121015 

II 

0.287 

-22.00 

Benetti et al. 

2014 

\ 

SSS120810 

Ic 

0.156 

-20.45 

Nicholl et al. i 

2014] 

I 

PSl-llap 

let 

0.524 

-20.54 

McCrum et al 

.1201 

3 

‘Silver’ sample: rest-frame g-band with bolometric correction 

SN2005ap 

Ic 

0.283 

-21.22 

Quimby et al. 

(2007 

7 

SCP06F6 

Ic 

1.189 

-21.56 

Barbary et al. 

200S 


PTF09cnd 

Ic 

0.258 

-21.34 

Quim 

by et al. 

(2011 


PTF09cwl 

Ic 

0.349 

-21.15 

Quim 

by et al. 

(2011 


PSl-lOky 

Ic 

0.956 

-21.24 

Chomiuk et a 






(2011) 



PSl-lObzj 

Ic 

0.650 

-20.32 

Lunnan et al. (2013) 

iPTF13ajg 

Ic 

0.740 

-21.50 

Vreeswijk et al. 


pH 1; 

*Pseudobolometric magnitude at maximum light; ^Described in 
the literature as a slowly-declining event 


2.1 Published PESSTO objects 


The first batch of SLSNe Ic classified by PESSTO were pre¬ 
sented and analysed by Nicholl et ah (2014). These were 


LSQ12dlf, SSS120810 and SN 2013dg. Each object exhib¬ 
ited spectral evolution typical of the class, despite their 
light curves being quite diverse. Another PESSTO SLSN, 


CSS121015, was studied by Benetti et al. (2014). This object 


was an extremely luminous SLSN II, but bore resemblance 
to SLSNe Ic in both the spectrum and overall light curve 


shape. Fitting by Nicholl et al. (2014), with magnetar- and 


CSM-powering, showed that for a given power source, the 
CSS121015 models occupied a similar region of parameter 
space to the SLSNe Ic. Including CSS121015 in our sample 
may help to clarify the existence of a link between normal 
SLSNe Ic, and some SLSNe II. 


2.2 LSQ14mo 

LSQ14mo was discovered rising steadily in LSQ observa¬ 
tions taken from 2014 Jan 12.2 UT. The transient is located 
at RA=10:22:41.53, Dec=-16:55:14.4 (J2000.0). A spectrum 
taken by PESSTO on 2014 Jan 31.2 UT was dominated by 
a blue continuum and O II absorption at around 4000 A, 
revealing it to be a SLSN Ic at a phase of ~1 week before 
peak luminosity. The spectrum was an excellent match to 
PTF09cnd (Quimby et al. 2011) at a redshift z ~ 0.25. A 


precise redshift of z = 0.253 was subsequently determined 
from narrow Mg II AA 2795,2802 absorption (Leloudas et al.| 


2014). PESSTO has collected extensive data on this target, 


which will be presented in full in a future publication (Chen 
et al., 2015, in prep.). 


2.3 LSQ14bdq 

LSQ14bdq was also discovered by LSQ during a long rising 
phase (>40 days), with the first confirmed detection occur¬ 
ring on 2014 April 5.1 UT, at a position RA=10:01:41.60, 
Dec=-12:22:13.4. A spectrum was obtained by PESSTO on 
2014 May 4.9 UT, showing it to be a SLSN Ic before maxi¬ 
mum light. The redshift was determined to be z = 0.347, 
initially by comparison with other pre-maximum SLSNe 
such as PTF12dam and PTF09cnd, and then 
through the detection of Mg II absorption 
2045). 


more precisely 
(Nicholl et al. 


2.4 SN 2013hx 


A hostless transient was first detected by SkyMapper 
on 2013 Dec 27 UT at coordinates RA=01:35:32.83, 
Dec=-57:57:50.6. It was given the survey designation 
SKYJ1353283-5757506. PESSTO observed the object on 
2014 Feb 20 UT after it had risen in luminosity for >30 
days. The spectrum showed Ha emission at z = 0.13, at 
which redshift the absolute magnitude was ~ —22, as well 
as broad features in the blue. It showed similarity to both 
SN 2010gx (a prototypical SLSN Ic; Pastorello et al.||2010) 
and CSS121015 (SLSN II; |Benetti et al.||2014| ). TheSN has 


been followed up by PESSTO and given the IAU name, 
SN 2013hx. A separate followup paper will present the full 
dataset (Inserra et al., 2015, in prep.) 


2.5 SLSNe from the literature 


The amount of data available for objects in the literature is 
highly variable. In some cases, they have only been observed 
in one or two filters; in others, they are at high redshift and 
the observed optical light corresponds to ultraviolet (UV) 
emission in the supernova rest-frame. High-redshift SNe also 
tend to have sparse spectral data, as they are fainter for 
observers. 

The objects in our sample have therefore been divided 
into two bins (‘Gold’ and ‘Silver’ samples), depending on 
whether they have good coverage at rest-frame optical wave¬ 
lengths. This can be seen in Table [l] The mean redshift for 
Gold objects is (z) = 0.22. All but two of the SLSNe at 
z < 0.3 have extensive photometry in observer-frame g,r, i 
and in most cases z filters, which at this redshift covers the 
rest-frame optical regime. This includes all of the PESSTO 


objects, the 5 low -2 

: SLSNe Ic from Inserra et al. (2013), the 

prototypical SN 2010gx ( 

Pastorello et al. 2010), 3 slowly de- 

dining SLSNe Ic (classified as SN2007bi-like; Gal-Yam et al. 

2009 Nicholl et al. 

2013; McCrum et al.|2014) and one fur- 

ther type II event (SN2008es Gezari et al. 2009). PSl-llap, 


at z = 0.524, falls in this group because of photometry in the 
NIR Pan-STARRSl y filter, which corresponds to rest-frame 
z-band at this redshift (McCrum et al.|2014 ). 


In the Silver sample (with a mean of (z) — 0.63), we 
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Figure 1. Complete set of rest-frame g-band (top) and griz pseudobolometric ( bottom ) light curves, after extinction and ^-corrections. 
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Table 2. Comparison sample 


Name 

Type 

Mgriz 

Reference 



Well-observed SNe in 

the literature 



SN1994I 

Ic 

-16.79 

Filippenko et al. (1995|), 

SN1998bw 

IcBL* 

-16.84 

Richmond et al. (1996 ) 

Patat et al. '(2001) 

SN1999ex 

Ic 

-16.84 

Stritzinger et al. (2002) 

SN2002ap 

IcBL 

-16.49 

Mazzali et al.|(2002 





Gal-Yam, Ofek &; S. 

hemmerj 




(2002) 



SN2003jd 

IcBL 

-18.19 

Valenti et al. (|2008a) 

SN2004aw 

Ic 

-17.24 

Taubenberger et al. 

(2006| 

SN2007gr 

Ic 

-16.36 

Valenti et al. |2008b|) 

SN2008D 

lb 

-16.24 

Sode^oeq^t al. (2008), 




Modjaz et al. 1 

[2009 

r 

SN2009jf 

Ic 

-17.34 

Valenti et al. (] 

2011) 


SN2010bh 

IcBL* 

-16.97 

Cano et al. J2C 

)U) 


SN2011bm 

Ic 

-17.63 

Valenti et al. ( 

2012) 


SN2012bz 

IcBL* 

-18.82 

Schulze et al. ( 

2014) 

\ 

SDSS II sample from 

Taddia et al. (2615;) 

SN2005hl 

lb 

-17.55 




SN2005hm 

lb 

-15.85 




SN2006fe 

Ic 

-17.04 




SN2006fo 

lb 

-17.26 




14475 

IcBL 

-19.46 




SN2006jo 

lb 

-18.25 




SN20061c 

lb 

-17.31 




SN2006nx 

IcBL 

-19.71 




SN2007ms 

Ic 

-16.95 




SN2007nc 

Ic 

-16.82 





* Associated with observed gamma-ray burst 


have two more objects from the Pan-STARRSl Medium 
Deep Survey - PSl-lOky at z ~ 0.9 ( Chomiuk et al.|[2011 ) 
and PSl-lObzj at z ~ 0.6 ( Lunnan et al.|20l3 ) - while most 
of the others featured in the original Quimby et al. (2011) 
sample that defined the SLSN Ic class. Although several of 
these objects are at redshifts 0.25 < z < 0.3, the reddest 
available photometry is in the R-band, which corresponds 
to rest-frame B/g. SCP06F6 (z = 1.189) has HST i and 
z photometry ( |Barbary et al. 2009), which corresponds to 
rest-frame emission between the u- and g-bands. The final 
Silver object is iPTF13ajg (Vreeswijk et al. 2014), which 
has excellent photometric and spectroscopic coverage, but 
at z = 0.74 this mostly probes rest-frame UV. This means 
that for these objects we must rely on an estimated correc¬ 
tion to obtain bolometric light curves. 


3 BOLOMETRIC LIGHT CURVES 

To analyse the light curves of our SLSNe in the most ho¬ 
mogeneous way possible, we constructed two sets of light 
curves: rest-frame g-band magnitudes, M g , and pseudobolo- 
metric light curves covering rest-frame SDSS griz filters. The 
first step was to apply R-corrections to the observed magni¬ 
tudes, to transform them to the rest-frames of our objects. 
These were determined as follows. We calculated synthetic 
photometry, using the lRAi0package CALCPHOT, on all avail¬ 


able spectra, for each filter in which the SN was observed. 
Spectra were then corrected to rest-frame using DOPCOR, in¬ 
cluding a correction to the flux per unit wavelength by a fac¬ 
tor of 1 + z, and new synthetic magnitudes were calculated. 
The R-correction at the epoch of a given spectrum is simply 
the difference between the rest-frame and observed synthetic 
magnitudes. These corrections were then linearly interpo¬ 
lated to the epochs with photometry. For most cases we 
simply corrected g 0 bs -A #rf etc., aside from the following: 
LSQ14bdq (r ohs -A ^ RF ); SN 2005ap, PTF09cnd, PTF09cwl 
(Robs -A ^rf); PSl-llap (z 0 bs -A #rf; R-corrections taken 
from |McCrum et~aT. 2014]); PSl-lObzj, iPTF13ajg (Gbs -A 
#rf); PSl-lOky (2 0 bs — Prf); SCP06F6 (z 0 bs,^obs -a #rf). 
For the objects of [Quim by et al.| (2011), sometimes only one 
spectrum was available - in this case we also used the spec¬ 
tra of SN 2010gx from Pastorello et al. ( |2010| ) (a spectro¬ 
scopically typical event, with good temporal coverage), after 
artificially placing them at the desired redshift. Magnitudes 
were also corrected for Milky Way extinction according to 


the dust maps of Schlafly Sz Finkbeiner (2011), though host 
reddening was assumed to be negligible. 

For the Gold sample, the bolometric light curve was 
then calculated using this rest-frame photometry from the 
g to z filters. For LSQ14mo, the z -band was estimated us¬ 
ing the colours of SN 2013dg, to which it has a very similar 
light curve in the grz-bands, while the z-band magnitudes 
for PSl-llap and LSQ14bdq were taken from PTF12dam, 
as in |McCrum et ah (2014). Most objects have good cov¬ 
erage before/around peak in one filter only. In these cases, 
the colours were assumed to be constant, with values from 
the first epoch that had multi-colour photometry available. 
This is a reasonable assumption, as the colours show lit¬ 
tle evolution until 1-2 weeks after maximum light ( |Inser ra 


et al. |2013). A spectral energy distribution (SED) was then 


constructed by converting these magnitudes into flux at the 
effective wavelength of each filter. For all objects, the flux 
was set to zero blue-wards of the g-band and red-wards of 
z. The luminosity, L griz , was then calculated by integrating 
this SED over wavelength, and correcting for distance to 
the SN assuming a cosmology Ho = 72kms _1 , Qm = 0.27, 
and Qa = 0.73. This follows the procedure of |Inserra et al.| 
(2013). We express this in terms of pseudobolometric mag¬ 


nitude, using: 


M ( 


3rlz ~ 2 ' 51Ogl0 l 3.055 xlO 35 


(1) 


based on the standard definition of bolometric luminos¬ 
ity. Our complete set of g-band and pseudobolometric light 
curves are shown in Fig. [l] 

However, a different process was needed to derive the 
M© light curves of the Silver objects. It was possible to find 
the average bolometric correction, M gr i Z —M g , as a function 
of time (we define t — 0 as the epoch of maximum luminos¬ 
ity), for the Gold sample, and apply this correction to our 
other objects. This correction is shown in Fig. [2] Unfortu¬ 
nately, there is significant scatter, but it is clear that the 
bolometric correction becomes more negative as a function 
of time. This is expected; as the SNe cool, bluer wavelengths 


3 Image Reduction and Analysis Facility (IRAF) is distributed by 
the National Optical Astronomy Observatory, which is operated 
by the Association of Universities for Research in Astronomy, Inc., 


under cooperative agreement with the National Science Founda¬ 
tion. 
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Figure 3. Interpolating the light curve of SN 2010gx. The times 
at which the dashed line intersects the polynomial fit give the 
exponential rise and decline times. 


Figure 2. Estimating the time-dependent bolometric correction 
for a typical SLSN. Our best fit is M gr i z —M g = 0.90 — 0.012£ 
for t > 0, though the gradient can vary from this by a factor ~2 
for individual objects. Uncertainty in the y-intercept term has no 
effect on our analysis. 

fade faster at late times. We fit only points where t > 0, 
assuming a constant correction before this. Our best fit is 
M g riz—M g — 0.90 — 0.012t for t > 0 and M gr i Z —M g — 0.90 
for t < 0. The uncertainty in this correction is ^0.5 magni¬ 
tudes by 50d after peak, hence at late phases the Mq light 
curves of Silver objects may become unreliable. However, 
note that only two Silver objects have data at this phase: 
PTF09cnd and PTF09cwl. Both of these objects closely re¬ 
semble PTF12dam in rest-frame g (but were not classified 
as being 2007bi-like), and this resemblance is preserved in 
the pseudobolometric light curves. Hence we feel justified in 
including these objects in our sample. 

4 LIGHT CURVE TIMESCALES 
4.1 Measurements 

Having constructed our pseudobolometric light curves, we 
proceeded to measure the rates at which these supernovae 
brighten to their maximum luminosity, and subsequently de¬ 
cline. We define: 

• T r i S e (the rise timescale): the time (t < 0) relative to 
maximum light (Lmax) at which L gr i Z — L max /e 

• Tdec (the decline timescale): the time (t > 0) relative to 
maximum light (L max ) at which L griz = L max /e. 

We make our measurements by fitting the light curves 
with low-order polynomials. Order four polynomials were 
found to give a good fit to all of our light curves for epochs 
t < 50 d. The fits were used to make a new estimate of 
the date when the pseudobolometric luminosity peaks. This 
tends to be later than the peak in y-band, as one would 
expect since the ejecta cool over time. We define the new 


peak, and then measure the quantities described above by 
interpolating the light curves with the polynomial fits. The 
method is demonstrated in Fig. [3] In some cases, the rise 
time had to be estimated by extrapolation using our poly¬ 
nomials. We consider this to be reliable for most objects, 
where we extrapolate only by a few days, but the rise time 
is poorly constrained for SNe 2007bi, 2005ap, and PSl-lOky. 
For slowly declining objects, the fourth order fit to the peak 
was not always a good fit at late epochs; for these objects, 
we made one fit to t < 50 d to estimate the peak and the 
rise time, as for the rest of our sample, and then measured 
the decline time by fitting another polynomial to only the 
post-maximum data points (fourth order and linear fits gave 
similar results). Our measurements are given in Table[3] The 
rise time of SSS120810 could not be constrained from the 
available data. 

4.2 Correlation 

We plot the rise times vs decline times for our sample in 
Fig. El The best-fit lines to our Gold and complete samples 
are calculated as follows: we represent each data point by a 
two-dimensional Gaussian, with a mean given by our mea¬ 
sured values of T r i se and Td ec , and standard deviations by the 
error bars. A Monte-Carlo method is then employed. A point 
is drawn at random from each probability distribution de¬ 
fined by these Gaussians, and we use standard Python rou¬ 
tines to calculate a straight-line fit, as well as Spearman and 
Pearson correlation coefficients, for the resulting rise-decline 
relation. This is repeated 10000 times. As can be seen in 
Fig. [5] the data are clearly correlated with high significance: 
Spearman’s rank correlation coefficient is 0.77±0.10 for the 
entire SLSN sample (0.75±0.11 for Gold sample only). Pear¬ 
son’s test gives 0.81 ±0.15 (0.84±0.14). The best-fit straight 
line to the data is Td ec = (1.65 ± 0.33)r r ise ± (7.38 ± 7.79) for 
the full sample, and Tdec = (1.96±0.46)r r i se ±(—0.10± 10.19) 
for the Gold objects. Although the formal best fit is different 
for the full sample compared to the Gold SLSNe only, the 
entire population is clearly consistent with a straight-line, 
and the two fits agree within the errors. It is no surprise 
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Figure 4. Rise vs decline timescales for SLSNe and normal stripped-envelope supernovae. The rise and decline times are clearly correlated, 
and in a similar way for both samples. The dashed black line gives the best linear fit to the entire SLSN sample (y = 1.65a; + 7.38), and 
the solid black line to the Gold sample only (y = 1.96a; — 0.10). The T14 Ic light curve was constructed by integrating the griz templates 
for SDSS SNe Ic from Taddia et al.| (2015|. 


that the gradient is greater than unity, as the light curves 
of other SN types (both Type la and core-collapse) rise to 
maximum more quickly than they decline. We also measure 
rise and decline timescales for the SNe Ibc in Table 13 They 
obey a similar correlation as for SLSNe, but with shorter 
timescales than their more luminous cousins. 


4.3 Models: overview 

To interpret our correlation, we use the synthetic light 
curve code described by Inserra et al. (2013). This code al¬ 
lows us to model the luminosity from a homologously ex¬ 
panding spherical ejecta with constant opacity and a cen¬ 
trally located power source, Li n (t). The light curve equa¬ 
tion, re-derived following the original Arnett (1982) paper 


but adapted for arbitrary Li n , is 


/•S.\(f) = <■ (t/Tm)2 f 2 Lin(i')— e^' /Tm ) 2 —, 

J o 1~rn 7~m 


( 2 ) 


where r m is the diffusion timescale (formally, the geomet¬ 
ric mean of the expansion and diffusion timescales; |Arne tt] 
1982). 

In the most basic case (for a fixed form of the power in¬ 
put term, L- in , e.g. an exponentially declining term for heat¬ 
ing by radioactive decay, or a central engine with a power- 
law decline), our model takes three parameters: a diffusion 
timescale, a power input timescale (which we call ri n ), and 
an overall energy scale (which affects the luminosity of the 
light curve, but not the shape). In general, the diffusion 
timescale is a function of ejecta mass, opacity and expan- 
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Figure 5. Fit parameters for the rise-decline timescale correlation. 


sion velocity: 


[2nMei V /2 
V P cv ) 


(3) 


where k is the opacity, M e j is the ejected mass, (3 ~ 13.8 
(for a wide range of plausible density profiles) is an inte¬ 
gration constant, c is the speed of light and v is a scaling 
velocity for homologous expansion (Ar nett|1980[ fl982 ). For 
a given opacity and velocity, the diffusion timescale thus al¬ 
lows us to derive the mass. Other authors have taken the 
observed rise times of SNe as an estimate of r m (most re¬ 
cently Wh eeler, Johnson &; Clocchiatti| 2014). This is a rea¬ 
sonable approximation for 5b Ni-powered light curves, where 
the decay time is well known, and is closely matched to 
the typical diffusion times (a coincidence that results in the 
high peak luminosities in Type la SNe). It is not surprising 
that the normal-luminosity hydrogen-poor sample shown in 
Fig. [4] obey a tight rise-decline correlation: the power input 
time is the same for ah these 56 Ni-powered SNe, and hence 
the diversity in both rise and decline times is driven by only 
one parameter: the diffusion time. The fact that the SLSNe 
obey such a similar correlation suggests that the diffusion 
time may also drive the correlation in these objects. 


However, for SLSNe the power input time, Tin, is un¬ 
known, and may span a wide range of values. If Ti n is very 
different from r m , it can have a large influence on the ob¬ 
served rise time, which is no longer a reliable proxy to 
T m . A better method here is to use the light curve width: 
we estimate that the diffusion time through the ejecta is 
Tm ~ (T r ise + Tdec) /2. This is explored in detail for the fol¬ 
lowing models, and will be important when we later attempt 
to estimate masses, in section [9] 


4.4 Models: 56 Ni and generalised exponential 
models 

As a first step towards investigating the SLSN parame¬ 
ter space, we generated an array of models with a hypo¬ 
thetical exponentially decaying power source (i.e., with the 
same functional form as 56 Ni decay, but for a variable life¬ 
time) . While this power source is not motivated by any pro¬ 
posed physical model, it aids in understanding the relevant 
timescales, by virtue of being the simplest possible scenario. 
This model takes only two parameters: the diffusion time 
(r m ) and the input time (ri n ), with L iri — Lo exp(— t/ri n ) 
(Lq is arbitrary). We varied the diffusion time between 10 
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Figure 6. The observed rise and decline timescales of our sample overlaid on a grid of diffusion models ( Arnett||1982| l, with a central 
heating term that decays exponentially with time. Input times of 10-20 d are needed to reproduce the SLSN correlation, which is then 
driven by variation in the diffusion timescale, as shown by the colour-bar. 56 Ni/ 56 Co models are also shown (grey curve). This roughly 
matches the normal-luminosity SNe Ibc, altho ugh t here is a slight offset, which depends on the efficiency of gamma-ray trapping (our 
model uses the trapping formalism of; Arnett] ( |l982| , with a velocity of 10000 kms -1 ). 


and 120 days, in steps of 10 days. A subset of these models 
are shown on Fig. [5] We see that, if SLSNe are powered by 
some universal, exponentially declining process, such a pro¬ 
cess must have a lifetime of ~ 10-20 days (the smaller points 
shown are spaced by 1 d in T in ). However, if the timescale 
for power input is variable, SLSNe in the top right can have 
longer timescales (shown by the red curve and larger points, 
with Ti n — 50 d). Most importantly, this figure shows that 
for this simple model, the correlation in rise and decline 
timescales is driven by r m , as seen in the colour scale. 

Figure [7] shows the relationship between diffusion time 
and measured rise/decline times for a full grid of mod¬ 
els, with T m ,Tin — 10-120 d, in steps of 10 days. Clearly, 
there is no straightforward way to deduce the diffusion 
time directly from either the rise or the decline. However, 
we find a strong correlation when the quantities are com¬ 


bined: T r i S e + Tdec ~ + Tin. This holds over the full range 

in Tm and Ti n . Hence if we know the timescale of our expo¬ 
nential power source, we can accurately recover the diffusion 
time by measuring the rise and decline timescales. Measure¬ 
ment of the ejected mass through light curve fitting has been 
applied for many years to SNe Ibc. Since the input timescale 
is known for 56 Ni-powered SNe, the light curve width is typ¬ 
ically a measurement of the diffusion time (e.g. Arnett [ 1982 


Drout et al.|2011| Valenti et al.|2008a ). For SLSNe, most of 


the best-fitting exponential models have Ti n ~ r m . Therefore, 
from our relation between the 4 important timescales, an es¬ 
timate of the diffusion time is given by T m ~ (T r i se + Td ec )/2. 


Figure [6] also shows the expected rise-decline curve for 
models powered by 56 Ni-decay. In this case, the only im¬ 
portant variable is the diffusion time ( 56 Ni mass sets the 
overall luminosity, but not the light curve width). We see 


Diffusion timescale, r m (days) 
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Figure 7. The relationships between various timescales in the 
simple exponential model. While neither the rise nor decline 
timescale is a good tracer of the diffusion time (top), the light 
curve width is determined by the sum of the power input and dif¬ 
fusion times, such that r m can be deduced from observed r r i se and 
t dec for a given Ti n (bottom). 


that this model does predict the correlation exhibited by 
SNe Ibc, which must be controlled by the diffusion time 
as mentioned in the previous subsection. Although some of 
fast rising and fast decaying SLSNe Ic lie close to this 56 Ni 
decay curve, that power source has already been ruled out 


for these. As discussed in Chomiuk et al. (2011), Quimby 


et al. (2011) and Pastorello et al. (2010), the peak luminos¬ 
ity means that the bb Ni mass would have to be greater than 
or similar to the total ejecta mass. Such expanding balls of 
56 Ni are unphysical and ruled out by the observed spectra. 


4.5 Models: magnetar 

In one of the most popular models, SLSNe are powered by 
a central engine which re-shocks the ejecta after it has ex¬ 
panded to large radius, thus overcoming adiabatic losses. In 
the magnetar spin-down model, the energy source is the ro¬ 
tational energy of a millisecond pulsar, which is tapped via a 
strong magnetic field. It is generally assumed to radiate with 


Figure 9. Same as Fig. [7] but for our magnetar model grid. The 
models shown have diffusion times of 10-100 d (steps of 10 d), 
B = {1,3,5,7,9} x 10 14 G, and P = {1,3,5,7,9} ms. Only models 
with peak luminosity greater than 3 X 10 43 ergs -1 are plotted. 


the functional form of a magnetic dipole (Kasen & Bildsten 


2010 ): 


-'’magnetar 


(t) = 


l 


p (1 + t/r v y 


erg s 


(4) 


where E p (the rotational energy, I ns Q 2 / 2) and r p are deter¬ 
mined by the spin period, P , and magnetic field, B, of the 
magnetar. The energy input timescale for magnetar spin- 
down is given by r p — 4.75B^ 2 P^ s d, where P 14 is the mag¬ 
netic field in units of 10 14 Gauss and P ms is the initial spin 
period in milliseconds. The shortest possible rotation period 
(corresponding to the largest energy reservoir) is P ~ 1 ms; 
any shorter and centrifugal forces woul d lead to b reakup. 
Galactic magnetars have B < 10 15 G (e.g. 


Davies et al. 


2009 


Table 3). This combination makes it difficult to achieve spin- 
down timescales r p < 0.1 d. 

We ran a grid of magnetar models, uniformly varying 
Tm (in steps of 10 days), B and P (in unit steps of 10 14 G and 
milliseconds, respectively), which we compare to the data in 
Fig. i However, it is not obvious that we can vary B and 
P independently. The high magnetic field is likely generated 
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Figure 8. We overlay a grid of magnetar-powered diffusion models with different spin period ( P ), magnetic field (B) and diffusion time 
(r m ). All are varied in uniform steps. Only models with L pea k > 3 x 10 43 ergs -1 are plotted. The top panel shows how increasing the 
diffusion timescale, r m , traces out the correlation we see in our data, and that r m is the parameter most strongly driving the diversity 
in our light curves. The lower panels show the effect of varying P and B. The bottom left region (where normal SNe Ic reside) is difficult 
to reach with magnetar models, as very high B is required. SLSN he along the lower right edge of the magnetar distribution, where P 
and B are least extreme (for given r m ). 
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Figure 10. The distribution of the power input timescales at 
maximum light, for the magnetar models in Fig. [9] We caution 
that this is for a uniform distribution in B and P, whereas we 
do not know what initial spins and magnetic fields magnetars are 
likely to form with. In particular, models with Ti n < 10 d require 
high magnetic field and fast spin. 


by a dynamo mechanism during core-collapse (Dunc an fc 
Thompson |1992 ), as a primordial B-field in the progenitor 


core would couple its angular momentum to the envelope, 
braking the core and likely precluding formation of a mil¬ 
lisecond pulsar at collapse. The simplest assumption for the 
dynamo mechanism is that a constant fraction of the ro¬ 
tation energy, E p oc P -2 , is converted to magnetic energy, 
Fmag oc B 2 . In this case, we would have B oc 1/P. Using this 
approach (with uniformly distributed P), rather than uni¬ 
formly distributed B and P, results in no significant effect 
on the distribution of rise and decline timescales in Fig. [5] 
Bearing in mind that we do not know the initial distribu¬ 
tion of spin periods (and hence B-fields), we will continue 
our analysis assuming that B and P are independent pa¬ 
rameters for simplicity. 

Only models that are brighter than the faintest SLSN 
in our sample (PTFlOhgi) are plotted. We produced many 
more models of course, but the only ones that are relevant 
are those producing luminosities of the same order as the 
SLSNe. The models shown display a correlation in rise and 
decline timescales similar to that observed in SLSNe, al¬ 
though a small number of models have very slow declines 
relative to the rise time. The colour map shows that increas¬ 
ing Tm drives the models to longer rise and decline times, 
tracing out the observed correlation. We also investigate the 
effect of B and P. The P-field influences the rise time, but 
has little effect on the decline. P~10 15 Gis needed to reach 
Trise < 10 d. This could explain why no SLSNe are seen with 
such short rise times. The corollary also holds: if we were 
to observe SLSNe with rise times less than about 10 days, 
it would preclude the magnetar model would struggle to ex¬ 
plain them. Models with long P tend to rise quickly, but this 
is not a necessary condition, unlike the constraint on B. 

Most of our SLSNe lie below the magnetar grid, but are 
consistent (within the errors) with lying along the locus of 
points on the sharp lower-right edge of the model distribu¬ 


tion. This locus corresponds to weaker B, and P ~ 2-5 ms. 
If SLSNe are powered by millisecond magnetars, the less- 
extreme magnetic field may account for why the data pref¬ 
erentially lie at the lower edge of the model distribution. 
Overall, magnetar models satisfactorily reproduce the basic 
trend we see in the data. For a sensible range of parameters, 
diffusion time dominates the diversity in timescales, just as 
we saw previously for the simplest models with a central 
(exponential) power source. 

Although most of the energy injection in this model 
takes place within a time r p , the power input timescale 
is not a constant, unlike in exponential models. Using the 
definition n n = \Li n (t)/(dL/dt)\, we find that for magne¬ 
tars, Tin(£) = \{t p + £), whereas for radioactive sources, 
Tin is simply the lifetime of the nucleus. In setting the 
peak width, the important timescale (in addition to r m ) is 
Tin(tpeak) = §(T p +£ pea k)- Using this definition, we recover al¬ 
most exactly the same correlation between the 4 timescales: 
Trise+ Tdec + Tin (tpeak). This is shown in Fig. [9] Does 

this mean that we can assume T m « (T r i se + Td ec )/2 in this 
case also? Fig. 


10 


shows the distribution of Ti n (t pe ak) for our 
magnetar grid. The mean is 19 days (standard deviation: 18 
d). Of all the models, 76% have Ti n between 5 and 30 days 
at peak; however those with Ti n (t pe ak) < 10 d need both fast 
rotation and strong magnetic field. This corresponds to the 
bottom left region in Fig. [8] where we do not have any ob¬ 
served SLSNe. Therefore the relevant models have timescales 
mostly in the range 10-30 d, with a tail extending to many 
tens of days. This is good news, as the assumption that 
Pn ~ T m is thus also reasonable for magnetar models. Hence 
we conclude that T m « (T r i se + Td ec )/2, for a range of sensible 
models with central power sources. 


4.6 Models: CSM interaction 


Alternatively, we can attempt to fit the observed correla¬ 
tions with parameterised models in which the ejecta collide 
with a dense CSM. This has been another popular model in¬ 
voked to explain SLSNe( Ginzburg &; Balberg||2012 Cheva¬ 


lier fe Irwin|2011 Moriya & Maeda 2012). As in Nicholl et al. 
(2014) (which followed the physical treatment of Chatzopou- 
los, Wheeler & Vinko 2012), we consider the limit of zero 


expansion velocity and luminosity input by strong external 
shocks. This model has many free parameters: ejected mass; 
CSM mass, radius and density (and density profile); explo¬ 
sion energy; 56 Ni mass. We ran a grid of models with fixed 
explosion energy (10 51 erg; this mostly just affects the peak 
luminosity) and inner CSM radius (10 12 cm; the light curve 
is quite insensitive to this parameter), and assume no 56 Ni. 
The CSM is taken to be a spherically symmetric shell of 
constant density - its mass and density therefore determine 
its radial extent. We varied the ejected mass (M e j) in unit 
steps from 1-25 M 0 , and CSM mass (Me sm) in similar steps 
from 1-15 M 0 , with the additional restriction Mcsm < M e j. 
The CSM density was initially set to pcsm = 10“ 12 gem -3 . 

The results are shown in Fig. EH At a given M c sm, 
increasing M e j increases both T r i se and Td ec , in the same 
way as the diffusion time, r m , did for our magnetar models. 
Illustrative lines of fixed Me sm have been marked - the ob¬ 
served rise-decline correlation is best reproduced by models 
with McsM>M ej /2. This was a common feature of models 
presented by Nicholl et al. (2014), and in fact is probably a 
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Figure 11. Same as Fig. [8] but for interaction-based models. In this case, we vary the ejected mass and CSM mass, and CSM density. 
We find that the shape of the rise-decline distribution is highly sen sitive to the CSM d ensity - the obse rved correlation is approxi mately 
recovered for pcSM = 10 -12 gcm -3 . All of the light curve fits by |Nicho ll et al.j {2014) (and many by jChatzopoulos et al. 2013) found 
PCSM ~ 10“ 12 gem -3 . Increasing the ejecta mass moves light curves along our correlation, while increasing the CSM mass primarily 
affects the rise time. It can be seen that the correlation is best reproduced by the subset of models with McsM~M e j/2 - this was a 
general property of the light curve fits by |Nicholl et ah'] (2014). Models with lower CSM mass rise too quickly for a given decline rate. If 
ejecta-CSM interaction does power all SLSNe, our rise-decline correlation puts tight constraints on the progenitor systems. 


requirement for SLSNe powered by interaction, as the CSM 
mass must be an appreciable fraction of the ejecta mass to 
thermalise the bulk of the expansion kinetic energy. 

More restrictive, but perhaps more interesting, is the 
effect of CSM density on our light curves. The rise-decline 
relation is similar to our data for pcsM = 10 -12 gcm -3 . In 
Fig. |12[ we show the effect of varying pcsM for a few repre¬ 
sentative models. Marked on the figure is the approximate 
slope of the observed correlation for SLSNe. Changing pcsM 
to 10 -11 gcm -3 or 10 -13 gcm -3 moves our models to re¬ 
gions of the rise-decline plot far from where our data reside. 

Why does this density have such a strong effect on the 
ratio of rise and decline times in our models? The dom¬ 
inant energy source is the forward shock from the interac¬ 


tion, which deposits heat as it propagates through the CSM. 
At some point, the shock breaks out of the CSM shell, and 
can no longer contribute energy (a similar effect occurs with 
the reverse shock in the ejecta, but the forward shock turns 
out to be dominant in most cases); at this point our light 
curve usually peaks. The time taken for the forward shock to 
propagate through the CSM decreases with increasing den¬ 
sity (Chatzopoulos, Wheeler & Vinko 2012 equation 15). 
The model subsequently declines, as the stored energy dif¬ 
fuses out of the CSM, on the characteristic CSM diffusion 
time. This timescale increases with the CSM density, so 
that models with earlier peaks fade more slowly. Thus, for 
all other parameters fixed, there is an inverse relationship 

































































Diversity and ejected mass in SLSNe Ic 15 




Figure 12. Rise vs decline timescales for synthetic light curves 
powered by ejecta-CSM interaction. The relation is strongly de¬ 
pendent on the density of the CSM. Models with pcsm ~ 
10 —12 gem -3 (shown in Fig. EH correspond best to our data. 
This can also be seen here in the lower panel, where we com¬ 
pare the SLSN LSQ12dlf to 3 models with M e j = Mcsm = 5 Mq, 
and varying CSM density. The light curves peak when the for¬ 
ward shock from the ejecta-CSM collision breaks out of the CSM, 
and the subsequent decline is controlled mainly by the diffusion 
time in the CSM. Denser CSM results in faster shock propagation 
(shorter rise) and slower diffusion (longer decline), giving the in¬ 
verse relationship between rise and decline times apparent in the 
top panel. 


between rise and decline timescales as we vary pcsm , as seen 
in Fig. [12] 

Clearly, if SLSNe are powered by interaction with a 
dense CSM, our observed rise-decline relationship can place 
narrow constraints on the range of CSM densities present. Of 
the 6 SLSN light curves fit with CSM models by |Nicholl et al.| 
(2014), 5 had densities in the range —12.54 < log 10 pcsm < 
— 11.74 (no convincing fit was found for the final object, SN 
201 Ike). It seems contrived that virtually all H-poor SLSNe 
would have such similar circumstellar densities, particularly 
when modelling indicates that a range of densities can gener¬ 


ate the observed peak magnitudes (e.g. Chatzopoulos et al. 


2013 who also fit H-rich events). Three possibilities exist. 


The first and most obvious is that ejecta-CSM interaction is 


not the power source in SLSNe Ic. Alternatively, our simple 
model may not be a good description of interacting SLSNe 
(for example, the shape of the CSM density profile may 
be important, and not a uniform shell). One very impor¬ 
tant weakness in this analysis is that the interaction mod¬ 
els of Chatzopoulos, Wheeler &; Vinko|J2012 ) and Nicholl 


et al. (2014) treat the shocks following | Chevalier &; Frans- 


(1994), whose derivation was for Mcsm <C M e j, and it is 


unclear how the picture changes for massive CSM. Finally, 
some process in the evolution of SLSN progenitors might 
somehow be capable of consistently producing circumstel¬ 
lar environments within this density range. The homogene¬ 
ity in the spectral properties of SLSNe would then result 
from the similar physical conditions in the CSM. The last 
of these possibilities is intriguing, and determining this pro¬ 
cess could prove an important clue to understanding what 
kinds of stars produce SLSNe. However, observations of SNe 
known to be powered by CSM interaction (SNe Iln) show 
huge diversity and variation in their observed characteristics 
and inferred physical configurations. In any case, CSM mod¬ 
els will have to be able explain our observed correlation, if 
we are to continue to accept them as valid model for SLSNe. 


5 GENERALISED LIGHT CURVES AND 
PEAK LUMINOSITY 

The correlation in rise and decline timescales, presented in 
the previous section, was found to be the same for SLSNe Ic 
and for normal-luminosity, hydrogen-poor core-collapse SNe. 
This suggested that the two populations have the same basic 
light curve shape (and we inferred that the slower evolution 
of SLSNe was due to longer diffusion timescales). To fur¬ 
ther investigate the relationship between SLSN and normal 
SN Ic light curves, we construct a generalised light curve 
for each type of SN. We do this simply by taking the area 
in magnitude-time parameter space that contains all of the 
light curves in each sample (SLSNe of Type II are excluded 


from this analysis). This is shown in Fig. 13 where we also 
include the SN Ibc light curve template from |Taddia et ah] 
(2015), and a number of SNe Ibc with unusually high lumi¬ 


nosity. As expected, the SLSN light curves are much brighter 
and broader than typical SNe Ibc. In the unsealed SN Ibc 
light curve, the decline rate starts to change at around 30 
days after peak. This is because the ejecta are becoming op¬ 
tically thin, and the luminosity begins to track the decay of 
56 Co. We do not typically see his behaviour in SLSNe, for 
two reasons. For most SLSNe, the decline after peak is too 
fast to be compatible with realistic 56 Ni-driven models with 
such a high peak luminosity, as discussed in Sect. |4.4| If they 
do contain some small amounts of 56 Co (comparable to that 
in SNe Ibc), this is masked by the bright luminosity source 
that powers the peak. The more slowly declining SLSNe (SN 
2007bi-like), on the other hand, do match the 56 Co decay 
rate. If they were radioactively powered, the long rise times 
and broad peaks would allow them to join smoothly onto 
the 56 Co decay tail after maximum. 

We find that a very simple transformation maps the SN 
Ic light curves onto those of the SLSNe: an increase along 
the y -axis by 3.5 magnitudes (a multiplicative factor of ~ 25 
in luminosity), and a broadening in time by a factor 3. The 
most obvious interpretation of this correspondence is that 
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Figure 13. Generalised light curves for SLSNe compared with 
lower luminosity SNe Ic (normal and broad-lined). The blue area 
represents normal SNe Ic in our well-observed literature sam¬ 
ple. and those of|T addia et al.| ( |2015| >. SN Ic light curves can be 
mapped onto the SLSNe by a 3.5 magnitude increase in bright¬ 
ness and a stretch along the time axis by a factor of 3. Some 
broad-lined SNe Ic he in the magnitude gap, but tend to have 
narrower light curves than SLSNe. SN 201 lbm (V alenti et al.| 
2012) shows a normal Ic spectral evolution but has a light curve 
width comparable to SLSNe. 


the two sets of light curves are determined by the same un¬ 
derlying physics: the rapid expansion of shocked gas with 
small initial radius, heated by some internal power supply. 
The broader light curves of the SLSNe are indicative of a 
longer diffusion timescale (higher mass and/or lower veloc¬ 
ity) compared to the SNe Ic. The higher peak luminosity 
tells us that some additional energy source is heating the 
ejecta, compared to ~ 0.1 — 1 M® of 56 Ni in SNe Ic (this 
could be, for example, a millisecond magnetar). Indeed, this 
is the established theory explaining the diversity within the 
SN Ic class, where higher ejecta mass and 56 Ni mass result 
in broader and brighter light curves, respectively. 

If SLSNe Ic were powered by interaction of SN Ic ejecta 
with dense (and H-poor) circumstellar material, the light 
curve physics would be quite different: a combination of for¬ 
ward and reverse shocks in the ejecta and CSM, with a 
strong dependence on the various density profiles. Indeed, 
for the massive, optically thick CSM needed to generate 
super-luminous peak magnitudes, we should not actually see 
direct emission from the supernova ejecta until well after 
maximum light ( Benetti et al.|[2014 ). It would therefore be 
surprising to recover such a trivial transformation between 
normal and super-luminous Ic light curves. Circumstellar 



usual light curve shape with a plateau lasting for hundreds 
of days. 

One interesting question is whether there is a continuum 
of peak luminosities between normal and super-luminous 
Type Ic SNe. Since the discovery and characterisation of 


SLSNe, an apparent gap has been recognised. Richardson 


et al. (2014) have compiled large samples of SNe Ibc to de¬ 


termine absolute magnitude distributions in the standard 
Johnson iJ-band. Their study attempted to correct for bias 
and derive volume limited absolute magnitude distributions, 
however the targets do not all have enough data to deter¬ 
mine bolometric luminosities ( M gr i Z at maximum). There 
are some typing inaccuracies in the Richardson et al. sam¬ 
ple (e.g. 2006oz as a lb, and 2005ap as a type II), however 
the normal and broad-lined Ibc population does appear to 
have an upper limit to their peak brightness of around —18 
to —19, with the SLSNe sitting 3 magnitudes brighter. 

We compare the brightness distributions of our SLSNe 
and other SNe in Fig. [14] The SLSNe Ic peak at M gr i Z = 
—20.72 d= 0.59 magnitudes, while the normal SNe Ibc tend 
to peak at M griz — —17.03 ±0.58. We note here that jlnser ra 
|fe Smartt (2014) derived peak absolute magnitudes for a 
sample of 16 SLSNe (with a large overlap with our sam¬ 
ple) in a synthetic bandpass centred on 400 nm. They found 
(M400) = —21.86 ± 0.35. This seems to be more uniform 
than M g riz , though it is worth bearing in mind that in some 
cases our derived M gr i Z has been estimated using assumed 
colours, introducing additional uncertainty, whereas |lnserra| 


& Smartt (2014) had observations fully covering their syn¬ 


thetic bandpass for all of their objects at maximum light. 

Selection effects and bias in our sample mean we cannot 
definitively say whether there is an excess of hydrogen-poor 
supernovae with peak absolute magnitudes brighter than 
M g riz ~ —20 (i.e. SLSNe), or whether such events are the 
bright tail of a continuous magnitude distribution. The plot 
in Fig. [l4] is not meant be be a representative luminosity 
function, as one would require accurate relative numbers 
in either a magnitude or volume limited survey. We know 
definitively that SLNSe are rare and their relative rate with 
respect to normal SNe Ibc is of order 1 SLSNe per 3000- 
10000 SNe Ibc ( Quimby et al.|[20l3 McCrum et al.|[20l5 ). 
The SLSNe are obviously over-represented in our sample, 
relative to their rates of occurrence in nature, as we have 
compiled similar numbers of SLSNe and comparison SNe 
Ibc. To prove statistically that SLSNe comprise a separate 
population of events from the brightest ‘normal’ SNe Ic 
would require careful consideration of the selection factors in 
surveys such as PTF, PS1 and LSQ±PESSTO to determine 
unbiased relative numbers of SLSNe, SNe Ibc, and SNe Ic- 
BL and construct meaningful luminosity functions (e.g. as 
done for normal SNe in Li et al.|201 1) 

Nevertheless we can make some comments on the lu¬ 
minosity differences observed if we assume that the SNe 
Ibc we have compiled are fairly representative of the gen¬ 
eral population of such stripped envelope SNe. None of our 
spectroscopically normal SNe Ibc are observed to peak at 
M g riz ~ —19. However, some broad-lined SNe Ic do have 
peak magnitudes spanning the gap between normal and 
super-luminous SNe. Broad-lined SNe Ic are often high 56 Ni- 
producers, and some have been shown to be associated 
with observed gamma-ray bursts (GRBs) (Woosley & Bloom 


2006). The large Ni mass makes them brighter than typical 


SNe Ic. Fig. 13 shows two GRB-SNe, SN 1998bw (Patat et al. 


2001) and SN 2012bz ( Schulze et al.|20l4 ). The brightest Ic- 


BL in the sample of Taddia et al. (2015) is also shown. The 
object, SN 2006nx, was discovered at redshift z = 0.137, but 
if there was an associated GRB, it was not seen. SN 2006nx 
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actually has a similar peak magnitude to SLSNe, suggesting 
that it may in fact be a member of that class. However, its 
light curve is quite narrow not only compared to SLSNe, 
but also for such a luminous 56 Ni-powered SN. Its unknown 
spectral evolution precludes a robust answer to the question 
of whether it is physically related to the SLSN population. 


It is interesting that some SNe Ic-BL/GRB-SNe do lie 
in the gap between SNe Ic and SLSNe, as some authors have 
suggested that SLSNe and GRBs are in fact related. |Lun-| 


nan et al. (2014) argued that the two types of explosion oc¬ 


cur in similar low-metallicity environments. Lelou das et al.| 
(2015) also claimed that their host galaxies are similar, but 
also that those of SLSNe are more intensely star-forming 
than those of GRBs, perhaps implying more massive pro¬ 
genitors for SLSNe. Magnetar/engine-powered models have 
been invoked to explain both the high luminosities in SLSNe 
and the relativistic jets in GRB-SNe( |Thompson, Chang &; 


Quataert 2004; Metzge r et al.|20TT ). However, the required 


magnetar parameters in GRB models are more extreme than 
for SLSNe, with spin periods ~lms and magnetic fields 10- 
100 times stronger, in order to drive a jet that punches 
through the stellar envelope ( Bucciantini et ah||2009 ). Gal- 


Yam (2012) notes that some broad-lined SNe Ic, such as SN 


2007D (Drout et al. 2011) and SN 2010ay (Sanders et al. 


2012), reached peak magnitudes close to those of SLSNe, 


and likely required an additional energy source on top of the 
inferred ~ 1 M® of 56 Ni. Perhaps bright SNe Ic-BL represent 
events where the central engine enhances the luminosity as 
well as the kinetic energy of the explosion, but not to the 
same degree as in SLSNe, where the additional power source 
overwhelms the input from 56 Ni decay. In this framework, 
there could thus be a continuum in luminosity between the 
various subclasses of SNe Ic, depending on the properties of 
the engine. 


If, however, there is a significant gap in peak bright¬ 
ness, with very few objects intermediate between normal 
and super-luminous SNe Ic, it may be difficult to explain for 
circumstellar interaction models. For example, Type II su¬ 
pernovae show a broad but continuous distribution in peak 
magnitudes: from SNe II-P through Type II-L’s to bright 
Type Iln and SLSNe Iln. Some authors have argued that 
this hierarchy is driven by varying degrees of circumstellar 
interaction (Chevalier & Fransson 1994; Richardson et al. 


2002 


Gal-Yam||2012 ). A gap between SNe/SLSNe Ic could 


be more indicative of a threshold process - for example, 
only progenitors above some critical mass can form magne- 
tars, or undergo a particular instability (e.g. pulsational pair 
instability; Woosley, Blinnikov Heg er||2007 ). In the mag¬ 
netar model, there may also exist an observational ‘desert’ 
between cases where the spin-down power is sufficient to 
drive a stable jet (GRB), and cases where the spin-down is 
weaker and slower, forming a wind nebula and enhancing the 
late-time luminosity instead (SLSN). If the magnetar spins 
down very quickly but the jet does not break the stellar en¬ 
velope, there may be neither a GRB nor an enhancement in 
optical brightness (B. Metzger, private communication.) In 
summary, peak magnitude distributions of SNe from large, 
homogeneous samples are needed, and will provide an impor¬ 
tant constraint on the possible relationships between SLSNe 
and SNe Ic. 


Table 3. Measured properties and derived masses 


SN 

Trise (d) a 

Tdec(d) b 

^5169 

M e j /M 0 




(kms -1 ) c 


Super-luminous SNe 

SN2007bi 

SN2008es* 

32.1 

18.3 

84.5 

38.0 

11900 

01 1 +34.3 
0 I.I — 21.7 

q 7 + 3.0 

2.1 

SN2010gx 

15.2 

29.1 

10900 

4 -411:1 

SN2011ke 

15.2 

25.7 

10200 

q q + 1-9 
* "*- 1.5 

SN2011kf 

17.5 

28.5 

9000 

q 7+2.O 
-1.5 

SN2012il 

14.1 

23.2 

9100 

2 4+ 1 - 3 

z -^-i.o 

SN2013dg 

17.8 

30.7 

14000 

« 0+3.8 
°-o_2.9 

SN2013hx 

20.7 

33.6 

6000 

0 4+I.8 
O.^-i 4 

LSQ12dlf 

20.1 

35.4 

13700 

8-i±l:S 

LSQ14mo 

17.4 

27.3 

10200 

0 q + 1-9 
^• y -1.5 

LSQ14bdq* 

31.8 

71.2 

- 

20 . 4+46 

PTFlOhgi 

21.6 

35.6 

4800 

3-0+1 

PTFllrks 

13.2 

22.3 

18100 

4-4+1 

PTF12dam 

CSS121015 

37.9 

20.4 

72.5 

37.8 

10500 

10000* 

27 n+ 19 - 6 
Z ' - U -14.3 

6 . 5+1 

SSS120810 

- 

30.2tt 

11200 

5 . 7+1 

PSl-llapt 

35.3 

87.9 

- 

OQ O+30.3 
zy - z -19.6 

SN2005apt 

>11 

28.8 

- 

3 . 0 +I 

SCP06F6t 

28.8 

39.8 

- 

O 1+3-1 

U.l_2 7 

PTF09cndt 

PTF09cwlt 

32.0 

34.8 

75.3 

60.6 


99 9+21.5 

zz - z _14.3 

17.5+° 2 8 

PSl-10kyt 

13.7 

32.5 


4 1 + 4.0 

U.l—2 7 

PSl-lObzj 

18.4 

37.3 

130(H) 

7.8+1 

iPTF13ajg 

27.4 

62.0 

9100 

14.0++ 

Other H-poor SNe 

SN1994I 

6.8 

9.9 

10100 

n p;+ u - 2 

u.O-o 2 

SN1998bw 

12.2 

20.1 

26600 

5 . 3+1 

SN1999ex 

12.1 

17.8 

9300 

1 6 +0-7 
± -°-0.6 

SN2002ap 

8.1 

19.0 

20400 

2 9+ 2,8 

1.9 

SN2003jd 

9.6 

17.5 

16400 

9 q + l -5 
Z -°-1.2 

SN2004aw 

13.5 

28.2 

12100 

4 1+3-4 
z 1-1_2.4 

SN2007gr 

10.0 

17.6 

8400 

1 9+O.8 

1-Z —0.6 

SN2008D 

14.6 

24.1 

8700 

2.5+1 

SN2009jf 

14.2 

23.8 

10100 

2.8+1 

SN2010bh 

8.6 

16.3 

35000 

4 2+ 3 - 0 

SN2011bm* 

24.0 

57.2 

- 

12 . 7 ++ 

SN2012bz 

11.7 

18.4 

23000 

4-0+1 

SN2005hl 

14.3 

29.1 

5450* 

2.0++ 

SN2005hm 

11.5 

27.5 

9470* 

2-8+1 

SN2006fe 

11.3 

27.7 

5000* 

1-5+1 

SN2006fo 

14.1 

32.7 

10500* 

4.4+1 

14475 

7.3 

14.1 

18700* 

1 6 +1-2 

1 -U — 0.9 

SN2006jo 

7.3 

11.3 

14400* 

1 n+0-5 
i - U -0.4 

SN20061c 

9.3 

19.0 

9100* 

1 4 +1-1 
1-^—0.8 

SN2006nx 

7.1 

16.5 

15400* 

1 . 7 +I 

SN2007ms 

16.1 

34.6 

11400* 

5.6+1 

SN2007nc 

10.6 

22.2 

12700* 

2-6+:i 


Masses derived using equation [H] with k = 0.1 cm 2 g 1 and r m = 

( T rise + r dec)/ 2 . Error bars correspond to estimates with r m = 
t rise (lower) and Tm = T & ec (upper). a Characteristic time for SN 


to rise to maximum light, defined using r r i se = t(T pe ak/ e ) for 
t < tpeak5 b characteristic time to fade, r^ ec = t(L pea ^/e) for t > 


tpeakl c Velocity from minimum of Fe II A5169 absorption, 20-30 
d after maximum light; * Assumes v ~ 6000 km s —1 , based on SN 
2013hx spectrum; ^Velocity from the literature, not necessarily 
Fe II; 1 Assumes v = 10000 km s —1 ; **For mass estimate, we take 

Trise = 7-dec /I-6 (see Fig. [ 4 } 
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Figure 14. The distribution of magnitudes at maximum light. 
The scatter is quite low, with a standard deviation of only 0.50 
magnitudes (0.58 if we include the two SLSNe II). The mean peak 
magnitude for normal SNe Ic is —17.03. A few broad-lined events 
have luminosity comparable to the fainter SLSNe. The data are 
binned in 0.5 mag intervals. 


6 TWO TYPES OF SLSN IC? 


For most SLSNe, the data unambiguously exclude 56 Ni and 
56 Co decay as the main power source around light curve 
maximum: the 56 Ni-mass needed to power the peak is of 
the order of the total ejected mass inferred from light curve 
fitting, and moreover exceeds the limiting 56 Ni-mass in¬ 
ferred from the late-time luminosity (e.g. Inserra et al.|2013 ). 
However, Gal-Yam (2012) proposed that the events with 


broader light curves, such as SN 2007bi, are radioactively 
powered, and likely exploded as pair-instability SNe. These 
would be fundamentally different from the other SLSNe, 
having a different explosion mechanism and power source, 
and would be characterised observationally by longer light 
curve timescales. We now look to see if we can distinguish 
two distinct classes in our data. In fact, Fig. [4] shows a rel¬ 
ative lack of SLSNe with T r i Se ~ 25-30 d and with Tdec ~ 40- 
60d. To make this clearer, we plot histograms of r r i S e, Td ec , 
and the g-band decline in magnitude 30 days after maximum 
light. This is a proxy for the decline timescale that is much 
simper to observationally measure, and particularly useful 
for parameterising the decline in high-redshift SLSNe, which 
may not have good rest-frame coverage in the redder bands. 
The data are shown in Fig. |15| 


We apply the Dip Test of|Hartigan &; Hartigan (1985) 


to each parameter. The dip statistic, D, measures the max¬ 
imum difference between the empirical distribution and a 
unimodal distribution (chosen so as to minimise D ). A larger 
value of D indicates that the data are not well described 
by a unimodal probability distribution function (PDF). We 
test for multi-modality using a bootstrapping method. We 
construct 5000 random sets of length n, where n is the 
size of our SLSN sample (n = 19 for our full sample), 
drawn from a uniform PDF, and calculate D for each set. 
The probability p-value for the null hypothesis to be cor¬ 
rect (i.e. that the data are unimodal) is then given by 
p = iV(T>sLSNe < Dboot)/5000 (this is the fraction of ran- 




Figure 15. Histograms showing the rise (top) and decline (mid¬ 
dle) timescales (binned in 5 day intervals), and the g-band decline 
(binned in 0.25 mag intervals) in 30 days after maximum (bottom). 
While the distributions show some indication of bimodality by 
eye, applying Hartigan’s Dip Test ( |Hartigan fc Hartigan||T985| ) 
shows that this is not statistically significant. 
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Table 4. Significance testing for bimodality 


Property 

Sample 

Dip Statistic, D p-value 

Rise time 

Full sample 

0.060 

0.801 


SLSNe Ic 

0.069 

0.666 

Decline time 

Full sample 

0.067 

0.605 


SLSNe Ic 

0.077 

0.439 

A£30 

Full sample 

0.069 

0.722 


SLSNe Ic 

0.078 

0.357 

Estimated mass 

Full sample 

0.044 

0.992 


SLSNe Ic 

0.050 

0.974 


‘Full sample’ includes the 3 events with hydrogen lines in their 
spectra. 


dom sets, drawn from a uniform PDF, which appear to be 
less unimodal than our data). 

The results of the test (D and p) are shown in Table 
[4] In no case do we find statistically significant evidence for 
bimodality. Therefore we have no confidence in rejecting the 
null hypothesis, and cannot confirm the existence of two sub¬ 
classes (fast/slow) of SLSNe. The similarity in overall light 
curves (up to some stretch factor) presented in sections 4|5 


may actually be evidence in favour of a single class, and if so, 
likely a single distribution in timescales. We note, however, 
that we do see lower p -values when we remove the hydrogen- 
rich events from the sample - but not low enough to exclude 
a unimodal distribution of light curve timescales. Many more 
objects will need to be observed before a stronger statement 
can be made. If the SLSN population does turn out to be 
significantly bimodal when a larger sample is constructed, it 
could indicate either two separate explosion mechanisms, or 
two progenitor channels for the same explosion mechanism. 


7 SPECTRAL EVOLUTION 


Quimby et al. (2011) first presented pre-maximum spectra 


of a sample of SLSNe from PTF; the high degree of simi¬ 
larity enabled the authors to determine that these objects 
together formed a new class of supernova. Since Past orello| 
et al. (2010) showed the spectroscopic evolution of SLSN 


2010gx into a more typical SN Ic, ah H-poor SLSNe have 
been seen to follow this path, including slowly-declining ob¬ 
jects (N icholl et al.| 2013). Here we compile high signal-to- 
noise spectra of 11 objects from our sample, to construct 
the complete spectral evolution, shown in Fig.[l6] The data 
for these are from the references in Table fll and are avail¬ 


able from WISeREP (Yaron & Gal-Yam 2012) or via the 


a 

20 ] 


PESSTO data release through the ESO archive (Smartt 
et al.|2014 ) 


The pre-max spectra are extremely blue, peaking at 
around 2500A, and dominated by high-ionisation lines - 
particularly O II at around 4000A, seen mainly in absorp¬ 
tion. Blackbody fits at this phase give a colour temperature 
T co i ~ 15000 K. There is a major change at around maximum 
light: as the ejecta cool, the O II lines disappear (oxygen, 
with an ionisation potential of O I of 13.6 eV, is mostly neu¬ 
tral at T < 15000 K), leaving the optical spectrum largely 
featureless. Within 10-20 days after the luminosity peaks, 
broad P-Cygni lines of singly-ionised metals (with lower first 
ionisation potentials, ~6-8eV) emerge, mainly Ca II H&K, 
Mg II, Fell blends (particularly between 5000-5500A) and Si 
II. By 30 days, these lines are strong, and the spectra closely 


resemble a conventional SN Ic at around maximum light 
(though even at this phase, SLSNe are often bluer than SNe 
Ic at peak). This similarity, and the slightly bluer colours, 
are highlighted in Fig. [lT| 

The next lines to appear are O I A7775 and the Ca II 
near-infrared triplet, along with the semi-forbidden Mg I] 
A4571 emission. This line becomes dominant over the Mg 
II/Fe II P-Cygni somewhere around 30-50 days after max¬ 
imum, and can be difficult to disentangle from the allowed 
transitions in the region. Strong forbidden emission lines do 
not appear in the spectra until over 100 days after maximum 
light (slowly declining objects such as PTF12dam are still 
not fully nebular at 200d). In a normal Ic, [O I] AA6300,6363 
and [Ca II] AA7291,7323 are already prominent at 60-100 
days after peak. These do eventually become the strongest 


lines in our very late spectra (PTF12dam at 500d; Chen 
et al. •inny The late appearance of nebular lines in SLSNe 


may indicate densities higher than those in normal SNe Ic. 
However, the nebular transition point depends most strongly 
on the opacity of the ejecta, which is in turn determined by 
the ionisation of the dominant elements (opacity drops by 
several orders of magnitude when these recombine). There¬ 
fore the higher temperatures in SLSN ejecta could account 
for the slower evolution towards the nebular phase. 

In Fig. [IS] we plot the colour temperatures measured 
from blackbody fits to our spectra (using the whole observed 
range in wavelength), as well as fits to the photometry of 
some well-observed low-z SLSNe Ic. The temperature evo¬ 
lution is markedly different from other Type Ic SNe, includ¬ 
ing the energetic SN 1998bw. SLSNe display (and maintain) 
much higher temperatures before and around peak light, 
and cool much more slowly afterwards. This requires sig¬ 
nificant heating over an extended period of time (several 
tens of days); therefore the abnormally high UV emission 
in SLSNe is a consequence of sustained energy input rather 
than a shock breakout phenomenon. This could plausibly 
be explained either by central engine models (as shown by, 
e.g. 


Inserra et al. 


2013 


Nicholl et al. 2013), or by CSM- 


interaction models, provided the CSM is sufficiently ex¬ 
tended that the forward shock continues to propagate and 
deposit heat for tens of days after the supernova becomes 
visible. The temperature behind a strong shock is approxi¬ 
mately T s m 3/16 /c^ 1 , where v s is the shock speed, ra 

is the average particle mass and ks is the Boltzmann con¬ 
stant (McKee & Draine 1991). For singly-ionised, oxygen- 
rich material, the shock speed needed to reach T ~ 15000 K 
is v s ~ 12000 km s -1 . In 10 days, such a shock travels 
~ 10 15 cm, which is approximately the blackbody radius of 
the photosphere in SLSNe - thus the temperature evolution 
does seem to be consistent with interaction models. 

However, there are several features of SLSN spec¬ 
tral evolution which may be difficult to reconcile with 
interaction-powered models. In particular 


• We do not see narrow lines from slow-moving material 
in any SLSN Ic (though we do in CSS 12105). 

• Moreover, the dense, massive CSM required to match 
the luminosity should be extremely optically thick, so the 
spectral lines seen at early times should be from the 
outermost material. These lines are puzzlingly broad (> 
10000 km s -1 ) if they are not supernova ejecta. 

• The spectral evolution is very homogeneous among 
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Figure 16. Typical spectroscopic evolution of SLSNe Ic. The grey curves are representative blackbody fits to determine the colour 
temperature (shown in Fig. |18| >. The main line-forming ions are also marked. 

























































Diversity and ejected mass in SLSNe Ic 21 



18000 


3000 


4000 


5000 6000 

Rest wavelength (Angstroms) 


7000 


Figure 17. SLSN LSQ12dlf at around one month after maximum 
light, compared to the normal SN Ic 19941 at maximum. The 
dominant lines in the two spectra are identical. LSQ12dlf shows 
stronger continuum emission in the blue, despite being older. A 
SYN++ fit (from Nicholl et al. 2014) identifies the strongest lines. 
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SLSNe. Supernovae which are known definitively to be in¬ 
teracting with their CSM (i.e. SNe Iln) exhibit quite diverse 
spectra, generally showing gas with a range of velocities. 

Thus, the spectral evolution is more easily explained 
without circumstellar interaction. The most natural inter¬ 
pretation is that the ejecta from SLSNe have a similar com¬ 
position to that of SNe Ic, but because of the much higher 
temperatures, due to heating from a power source equivalent 
to several solar masses of nickel, we do not see normal SN 
Ic spectral lines until later in the evolution, instead seeing a 
very blue continuum and high ionisation lines around peak. 
The slow spectral evolution may be exacerbated by higher 
densities in SLSNe. 


Figure 18. The evolution of the colour temperature for SLSNe 
Ic and normal SNe Ic. Colours are the same as in Fig. |16| Tem¬ 
peratures are derived from blackbody fits to spectra (points) and 
photometry (lines). Photometric measurements tend to be hot¬ 
ter, as there is a lack of P-band data to constrain the fit at 
~ 3000 A, which is not a problem for the spectral fits. Crosses 
indicate UV spectra (A < 2000 A). The SNe Ic are 19941 (solid), 
2007gr (dotted) and 1998bw (dot-dashed). SLSNe appear to have 
a very slow colour temperature evolution prior to maximum light, 
with typical T co i ~ 15000 K. Around maximum light, the colour 
temperature drops rapidly, reaching ~ 10000 K around 10 days 
later. This is partly due to decreasing gas temperature, but the 
evolution is accelerated at this phase due to the emergence Mg II, 
Ca II and Fe II lines, which absorb flux in the blue. The SLSNe 
then cool at a constant rate for 50-60 days. Normal SNe Ic do 
most of their cooling soon after explosion, and reach a constant 
temperature (~ 5000 K) 10-20 days after maximum light. 


8 VELOCITY MEASUREMENTS 


The analysis in section ^ demonstrated that diffusion time 
in the ejecta is the most important factor in generating the 
observed diversity in SLSN light curves. As shown by equa¬ 
tion [3] this timescale is a function of ejected mass, velocity 
and opacity. Since all of our super-luminous objects (at least 
the hydrogen-deficient ones) have very similar spectral evo¬ 
lution, we expect that they have similar compositions and 
ionisation states, and hence similar opacity. This leaves us 
with a degeneracy in M e j and the expansion velocity, v. Ve¬ 
locities can be estimated using the spectra, allowing us to 
break this degeneracy, and to determine whether it is varia¬ 
tions in M e j or v that are most important between the SNe 
(super-luminous and normal) in our sample. We measure 
velocities using the absorption minimum of the Fe II A5169 
P-Cygni profile. This line profile is shown at t — 20-30 d in 
Fig. [19] 

For a few objects, we have sufficient temporal coverage 
and signal-to-noise to measure this over a period of more 
than 30 days, and see how the velocity evolves. This is shown 
in Fig. |20| The velocities are remarkably constant in time, 
declining by at most 2000km s -1 in the first 30 days after 
maximum light (and in many cases showing no clear decline 
at ah). This is quite discrepant with the velocity evolution of 
other SNe Ic, such as the sample from Valenti et al. (2012) 


spherical shock; Arnett 1980). Thus, as the ejecta expand 


shown in the figure. In normal SNe Ic, the ejecta expand in 
a roughly homologous fashion (a natural consequence of a 


and the photosphere recedes to deeper layers, we see more 
of the relatively slow-moving inner material, and hence mea¬ 
sure lower velocities. Constant velocities, on the other hand, 
are predicted by the magnetar models of Kasen V Bildsten 
(2010), who showed that for magnetars with rotational en¬ 


ergy greater than the kinetic energy of the supernova ejecta 
(i.e. > 10 51 erg), essentially ah of the ejecta is swept up into 
a dense shell of uniform velocity. This is an unavoidable con¬ 
sequence of the central overpressure, and thus far the clear¬ 
est observational test for engine-powered supernovae. The 
fairly hat velocity curves of our SLSNe are therefore con¬ 


sistent with the Kasen & Bildsten (2010) models. This was 


also pointed out by Chomiuk et al. (2011), in their study of 
SLSNe PSl-lOky and PSl-lOawh. Moreover, because of this 
slow evolution, we can use velocity measurements at 20-30 
days after maximum light as a reasonable proxy to the ve¬ 
locity at day 0, which is useful as for most objects we do not 
have good detections of Fe II until this phase. 

The velocity measurements from the Fe II line are given 
in Tab le [3| as well as velocities for the SDSS SN Ibc sample 
from Taddia et al. (2015), and the distributions for SLSNe 


and SNe Ibc are shown in Fig. [21] The velocity distribu¬ 
tions of SLSNe and normal SNe Ibc are almost indistin¬ 
guishable within the errors: the median velocity for SLSNe 
is 10500 kms -1 , with a standard deviation of 3100 kms -1 , 
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Figure 19. The Fe II A5169 P-Cygni profile in SLSNe at -20- 
30 days after maximum light. The minimum of the absorption 
trough gives an indication of the photospheric velocity. The line 
is clearly weaker (or absent entirely) in our SLSNe Type II (2008es 
and CSS 121015) - it does not become prominent in these objects 
until > 40 d after maximum light (Benetti et al. 2014j. 



Velocity (km s 1 ) 


while the median for normal SNe Ibc is 9800 kms -1 , with 
a standard deviation of 2500 kms -1 . The broad-lined SNe 
Ic, on the other hand, all have velocities greater than 15000 
kms -1 . The similarity between typical photospheric veloc¬ 
ities at peak light for super-luminous and normal SNe Ic 
indicates that the broader light curves and slower spectral 
evolution in SLSNe are not caused by a slower expansion. 
However, there is a possible caveat to this: since SLSNe 
have higher temperatures, and therefore higher ionisation, 
the photosphere may be formed further out in mass coordi¬ 
nate compared to normal SNe Ic, and therefore amidst faster 
moving ejecta. This would give a high photospheric velocity 
even if the bulk expansion was slower. Inserra et al.| (2013) 
and Nicholl et al. (2013) looked at the widths of [Ca II] and 



Figure 20. Velocity evolution of well-observed SLSNe Ic. All 
velocities have been measured from Fe IIA5169 absorption min¬ 
ima. The SN Ic velocities are taken from Valenti et al. (2012), 
who used the same method; these have been shifted downwards 
by 4000km s —1 for clarity of presentation. The velocity curves for 
SLSNe are seen to be much flatter than those of SNe Ic, where 
we see a rapid decline after maximum light. 


Mg I] emission lines in late spectra of SLSNe, which give an 
independent trace of the expansion velocity, finding typical 
values of 10000 kms -1 , fully consistent with our estimates 
from the Fe II P-Cygni lines. This gives us confidence that 
SLSNe expand at least as fast as other SNe Ic. 

Two of our objects, LSQ12dlf and SSS120810, do show 
significant velocity evolution. We also note that |Inserra et al.| 
(2013) found evidence of decreasing photospheric velocity in 
some SLSNe Ic, compatible with the predictions of simple 
spherical models, rather than dense shells. They presented 
Fe II velocities close to maximum light for a number of ob¬ 
jects, whereas here we had difficulty to reliably determine 
the Fe II profiles at such early epochs, when the line is very 
weak. The apparent velocity decline shown in their study 
also relied on measurements of different lines at different 
epochs, including emission line widths at late times. In fact, 
their Fe II velocities after 15-20 days from peak are flat¬ 
ter than at early times, and in reasonable agreement with 
our measurements here. However, if these objects do show 
an initial velocity gradient, it may indicate that the pro¬ 
cess causing the flat velocity evolution (for example, a sec¬ 
ond shock) has not yet terminated, or may suggest a tail in 
the ejecta distribution extending to high velocity. Another 
possibility is that the overpressure from a central engine is 


anisotropic (stronger along the polar axes; Bucciantini et al. 


2007). An observer might initially see high velocities from 


the faster polar region, with lower r m , before later seeing the 
contribution from the slower, higher-r m equatorial material. 
In any case, if peak velocities estimated from our measure¬ 
ments at 20-30 days are really lower limits, higher velocities 
at peak in SLSNe compared to normal SNe Ic would only 
serve to strengthen our conclusion: the broad light curves of 
SLSNe cannot be explained by low velocities. 
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Figure 21. Photospheric velocities measured from the lines 
shown in Fig. |19| compared with the velocities of SNe Ic 
at maximum light (the same objects in Fig. [ 4 J. Bin size is 
2000 km s —1 . Most SLSNe have velocities close to the mean value, 
« 10800 km s _ 1 , with a standard deviation of « 3000 km s _ 1 . The 
outliers are PTFllrks (fast) and PTFlOhgi (slow). The other slow 
object is CSS121015, though the Fe II line is very weak at this 
epoch (and not robustly detected in SN 2008es or SN 2013hx, 
the other SLSNe II). The median velocities are 10500 km s _1 , 
9800kms _1 and 20400 km s -1 for SLSNe, normal SNe Ibc and 
broad-lined SNe Ic, respectively. 


9 MASS ESTIMATES 

Having shown that velocity is not the parameter driving the 
diversity in SLSN evolution timescales, and if the opacity is 
similar for all of our objects, we are left with the ejected mass 
as the most important factor. Equation [3] can be rearranged 
to give 


M e j = 7.7 x 10“ 


( « ) 

1 v 1 

( Tm \ 

\0.1 cm 2 g -1 J 

1 kms- 1 1 

V days J 


(5) 


where all variables are as defined in section 4] and our veloc¬ 
ity measurements are taken from section [8 There is an im¬ 
portant caveat to our velocity measurements: we have mea¬ 
sured the velocity of the photosphere (I’phot), whereas the v 
in equations |3] and [5] are the ejecta ‘scale velocity’ (A rnett 
1980, 1982| ). There is no simple way to measure this scaling 
velocity, and it is unclear whether the photospheric velocity 
is a good representation of it. However, we will proceed with 
the assumption that v ~ t’phot- 

We take k = 0.1cm 2 g _1 as a fiducial value 


Inserra 


et al. (2013) showed that this is a reasonable approximation 


for opacity dominated by electron scattering, if the temper¬ 
ature behind the photosphere is ~ 10 5 K, since abundant 
species such as oxygen, carbon and iron will be roughly half- 
ionised. For fully-ionised gas, k = 0.2cm 2 g _1 . Shortly af¬ 
ter explosion, we expect full ionisation, but this drops as 
the expanding SN cools. The opacity will therefore vary 
over the duration of the light curve, but this is not taken 
into account in the Arnett (1982) formalism. Most au¬ 


thors use constant (essentially, time-averaged) opacities of 
ft = 0.07 — 0.1cm 2 g _1 for modelling normal SNe Ic. The 
ionisation fraction, and hence the electron-scattering opac¬ 
ity, may be higher for longer in SLSNe, due to the additional 


energy source, but this is unlikely to be by more than a small 
factor, simply because progressively more energy is required 
to remove successive electrons from an ion. From equation 
[5] increasing the opacity decreases the derived mass by the 
same factor. The lower limit on ejecta mass is found by as¬ 
suming full ionisation - this is therefore half of the mass 
estimated using k = 0.1 cm 2 g -1 . 

The final parameter to estimate is the diffusion 
timescale, r m . This was discussed in detail, both in general 
and for specific models, in section [4] but we recap here for 
convenience. It is common practice to use the rise time as 
an estimate of r m (for a recent example, see Wheeler, John- 


son & Clocchiatti| 2014 

). This i 

original formulation of 

Arnett 


but that derivation 
was only for 56 Ni-powered SNe ( &b Ni having an exponen¬ 
tial lifetime of 8.8 days, comparable to the diffusion time 
in SNe la). As the decay time of 56 Ni is fixed, the power 
input is the same for all normal Type I SNe, and any vari¬ 
ation in light curve timescales depends only on the diffu¬ 
sion time. However, SLSNe may have a range of power in¬ 
put times, for example magnetars with different spin-down 
times, which may be much longer or shorter than the dif¬ 
fusion time. While it remains true that the diffusion time 
should generally be represented in the shape of the light 
curve around peak, we lack a straightforward mapping be¬ 
tween Tm and either r r i se or Tdec- As a best estimate, we 
take Tm ~ (TVise +Tdec)/2, which is true for many scenarios, 
and use the rise and decline timescales as upper and lower 
limits on r m . The uncertainty in mass owing to the choice 
of timescale is a factor of < 2. 

More precise determinations can only be done through 
detailed modelling. Our goal here is not to derive the most 
exact masses, but to provide estimates for our sample in 
the most general and homegeneous way possible, thus ex¬ 
posing any underlying trends. The masses derived using our 
method are given in Table [3] One caveat we should add is 
that we have referred to the ejected mass, but this assumes 
that the entire diffusion mass associated with the light curve 
is supernova material; if SLSNe Ic are powered by CSM- 
interaction, the diffusion mass we are probing is a combina¬ 
tion of ejecta and CSM. However, most of our analysis in 
this paper indicates that SLSNe are governed by generally 
similar physics to SNe Ibc, i.e. rapidly expanding material 
being heated from the inside. We therefore propose that our 
derived masses are likely representative of the ejecta. 

The average ejecta mass in the SLSNe is 10.1 M© (stan¬ 
dard deviation: 9.0 M 0 ), compared to 3.1 M© (2.9 M©) and 
3.1 M© (1.3 M©) for the normal and broad-lined SNe Ibc, 
respectively. The median masses are 6.0M©, 2.5M©, and 
2.9 M© for the three samples. While we are sensitive to 
small-number statistics here, the peak of the broad-lined Ic 
distribution appears to be at higher mass than the normal 
SNe Ic, especially if the extreme outlier SN 201 lbm is ne¬ 
glected (giving a mean/median SN Ibc mass of 2.4/2.3M©). 
Taking only the SNe with an observed GRB counterpart 
(SNe 1998bw, 2010bh and 2012bz), the mean of 4.5 M© is 
intermediate between SNe Ibc and SLSNe. Fig. |22| shows the 
ejecta mass distributions for SLSNe, SNe Ibc and SNe Ic-BL, 
with the ejected mass higher for SLSNe by a factor of 2-3. Al¬ 
though the broad-lined SNe Ic appear to be skewed towards 
higher mass than normal SNe Ibc, overall they seem to eject 
significantly less mass than the SLSNe. This supports the 
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view of Leloudas et al. (2015), who suggested that SLSNe 


have more massive progenitors than LGRBs, but may be in 
tension with the results of Lunna rTet al.| |2015| ), who found 
that SLSNe do not trace host galaxy UV light (star for¬ 
mation) as closely as do LGRBs, implying older/lower-mass 
progenitors. 

There are clearly several very massive (M e j ~ 20-30 M©) 
SLSNe compared to the normal/broad-lined Ibc sample, 
for which only SN 201 lbm hints at a high-mass tail. Five 
SLSNe eject over 20 M©, yet there are only two objects in 
the 10-20 M© regime. To check whether this high-mass tail 
could be fully explained as a consequence of greater ioni¬ 
sation in SLSNe (e.g. due to the hard radiation held from 
a magnetar; Metzger &; Piro 2014), we rescale the masses 
to k, = 0.2 cm 2 g -1 in the bottom panel of the figure. In 
this case, the bulk of the SLSN sample have masses consis¬ 
tent with the normal SN Ic population. However, there is 
still a clear excess of events with M e j > 10 M©. Therefore 
even if many of the SLSN light curves can be explained as a 
consequence of high ionisation, there remains a substantial 
number of events that must eject significantly more mass 
than normal stripped SNe. 

The interesting question, then, is whether this indicates 
a separate population of high-mass SLSNe, arising from a 
different progenitor, explosion mechanism, or power source. 
Effectively, this is the physical interpretation of the simple 
observational result illustrated in Fig. [4] (and discussed in 
Section^ that the rise- and decline-time correlation visu¬ 
ally picks out two groups. We apply Hartigan’s Dip Test, as 
described in section [6] to our distribution of SLSN masses. 
We find that D = 0.050 and p = 0.974, showing that the 
most massive SLSNe are fully consistent with being the tail 
of a continuous distribution. However, many more objects 
will be required in order to test this more robustly. There 
may also be an observational bias here: broad light curves 
are easier to detect than fast ones, which could lead to the 
more massive objects being over-represented. 

The very large ejecta masses inferred for a few of our ob¬ 


jects may be surprising. For example, Nicholl et al. (2013) 


fit the light curve of PTF12dam with a magnetar engine 
in 16M© ejecta (for the same opacity used here). Based 
solely on the light curve timescales and assuming no par¬ 
ticular form for the power input, we have here estimated 
around 25 M©. Fortunately, this discrepancy can be easily 


explained. As described byTnserra et al. (2013), this magne¬ 


tar model assumes a kinetic energy FT ss 10 5r + 0.5(E mag — 
E r &d) erg, where E mag is the total energy input by the mag¬ 
netar and E ra d is that lost due to the radiation emitted from 
the SN. For SLSNe that evolve quickly, the magnetar spins 
down rapidly, and overwhelms the original kinetic energy, so 
the derived parameters are not very sensitive to the initial 
choice of 10 51 erg. However, for slow objects like PTF12dam, 
this is not necessarily the case (the fit by Nicholl et al. 2013, 
had FT = 1.69 X 10 51 erg), so the initial supernova explosion 
energy is important. In fact, for their kinetic energy and 
mass, the expected velocity is lower than observations by a 
factor of two. From this we draw several conclusions: 


• The explosion energy in PTF12dam exceeded 10 51 erg. 

• Velocity information can be an important constraint 
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Figure 22. The ejecta mass distribution for SLSNe and other 
H-poor SNe. SLSNe seem to arise from explosions ejecting > 2 
times the mass, on average, ejected by normal and broad-lined 
SNe Ibc, however, the mass distributions actually appear quite 
similar if the opacity in SLSNe is a factor 2 higher than in other 
SNe Ibc. A magnetar wind is a good candidate to increase the 
ionisation (hence opacity) in the ejecta ( |Metzger & Piro|20 14]> 


would be the scaling velocity of |Arnett]|1982[ though it is 
unclear how to derive this quantity from the spectrum). 

• Again using equation 3, for a fitted diffusion time, dou¬ 
bling the velocity means doubling the mass, so the same 
magnetar model with the observed expansion velocity (~ 
10000 km s -1 ) would have ejected twice as much material, 
consistent with our estimate here. 


We also note that Kasen & Bildsten (2010) repro¬ 
duced the light curve of SN 2007bi with a magnetar model 
and 20 M© ejecta, using an opacity k = 0.2cm 2 g _1 . In 
our estimates, we have used k = 0.1cm 2 g _1 . For a fixed 
explosion energy and diffusion time, halving the opacity 
means increasing the ejecta mass by a factor ~ 1.6 (Nicholl 
et al. 2014). Our estimate of 31.1 M© for SN 2007bi is thus 
quite consistent with the hydrodynamical magnetar model 


of Kasen V Bildsten (2010). This gives us confidence in our 


on the magnetar models of Inserra et al. (2013) (ideally this 


simple mass derivations. 

The masses also place constraints on explosion mod¬ 
els. A core-collapse explosion ejecting 20 M© might be ex- 
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pected to produce a black hole rather than a neutron star 
( Heger et al.| 2003). However, recent theoretical results have 
shown that the final fate of a massive stellar core (neutron 
star or black hole) is unlikely to be monotonic in mass (e.g. 
O’Connor & Qtt||2011| |Clausen, Piro &, Ott||2C)T5| ). In par¬ 
ticular, Dessart, O’Connor &; Qtt| ( |2012 ) have shown that 
in rapidly-rotating stars (potential GRB and SLSN progeni¬ 
tors), magnetars are produced more naturally. Nevertheless, 
if a black hole is formed in the core collapse, the central 
engine mechanism could still apply, in this case in the form 


of fall-back accretion (Dexter & Kasen 2013) rather than 


magnetorotational powering. It should be noted that the 


fall-back models considered by Dexter & Kasen (2013) only 


produced super-luminous light curves if the progenitors re¬ 
tained their hydrogen envelopes, in order to delay fall-back 
to later times. The progenitors of SLSNe are expected to be 


stripped of their hydrogen, though Nicholl et al. (2015) and 


Piro (2015) recently found evidence that they may exhibit 


an inflated He envelope at the time of explosion. 

Taking our upper limits on ejecta mass (from the decline 
timescales), the most massive objects may eject as much as 
40-60 M 0 , suggesting that pair-instability explosions should 
be considered as viable models. However quantitative com¬ 
parisons do not lead to comfortable agreement between the 
observed timescales and detailed models. Firstly, as shown 


by Gal-Yam et al. (2009), Kasen, Woosley & Heger (2011), 


Dessart et al. (2012) and Nicholl et al. (2013), H-poor PISN 
models with M e j> 100 Mq are needed in order to match 
the observed peak luminosities of SLSNe. This is a factor 
2 higher than our most optimistic mass estimates. There 
is also a larger problem, which is independent of our mass 
estimates. The 120 M© model of |Kasen, Woosley &; Heger| 
(2011) has Trise — 72 d and Td ec = 99 d, defined in the same 


way as for our objects (in these PISN models, the decline 
timescale is mainly set by 56 Co decay, which has a timescale 
of 111 d). A few of our objects have Td ec > 80d, but all have 
T r i S e shorter than the PISN model by at least a factor of 2. 


This is the same problem found by Nicholl et al. (2013) and 


McCrum et al. (2014) when investigating PISN models for 
PTF12dam and PS 1-1 lap. Our measurements here support 
their conclusion: that even the 2007bi-like SLSNe do not 
quantitatively match the pair-instability explosion models. 


broader and brighter. Broad-lined SNe Ic, including GRB- 
SNe, may span the magnitude gap between spectroscopically 
normal SNe Ic and SLSNe. 

• SLSNe Ic span a wide range of light curve timescales, 
but there is not yet statistically significant evidence for sepa¬ 
rate populations of fast- and slowly-evolving objects. It may 
be one continuous distribution. 

• The spectroscopic evolution requires sustained heating 
around peak, to maintain a temperature T co \ ~ 15000 K de¬ 
spite weeks of expansion. At around one month after maxi¬ 
mum light, SLSNe have cooled to temperatures comparable 
to normal SNe Ic at maximum, at which point the spec¬ 
tra show very similar photospheric-phase lines. This implies 
similar ejecta composition. 

• The temperature evolution is consistent with both 
central-engine and circumstellar-interaction light curve 
models. 

• The broad lines at all phases, lack of narrow lines in any 
object, similarity to SNe Ic, and the overall homogeneity of 
SLSNe spectra all argue against significant modification of 
the spectra by circumstellar material. 

• SLSNe and normal SNe Ic have similar photospheric 
velocities around maximum light. After 0-10 days from peak, 
this velocity stays remarkably constant in time for many 
SLSNe. This could be explained by a re-shock from their 
central power source, sweeping the ejecta up into a uniform 
shell ( Kasen &; Bildste"n||2()10 ). 

• For a given opacity, the ejected mass, derived from mea¬ 
sured Fe II A5169 velocities and simple estimates of the dif¬ 
fusion timescale, are on average > 2 times higher than the 
masses in normal and energetic SNe Ibc (although there is 
overlap). However, an alternative explanation is that the 
opacity in SLSNe is twice as high as in other SNe Ibc. Such 
an effect could possibly arise from a magnetar ionisation 
wind ( Metzger k, Piro| 2014). There is no statistical evidence 
for a bimodal mass distribution in SLSNe. 

• SLSNe may eject as much as several tens of solar masses 
of material. Our estimates are consistent with masses de¬ 
rived from magnetar light curve fits. The masses seem to 
be too low (and the light curve evolution too fast) to be 
consistent with pair-instability models. 


10 SUMMARY OF RESULTS 

For convenience, we here summarise our main findings: 

• SLSNe Ic typically have griz pseudobolometric magni¬ 
tudes in the range — 20 > M gr i Z > —21.5, with a mean of 
-20.72 and a standard deviation of 0.55. 

• Their light curves obey a fairly tight relationship be¬ 
tween rise and decline rates around maximum light, with 
Tdec ~ 2 T r ise- This is naturally produced by simple diffusion 
models, but the lack of scatter may be difficult to explain 
with models requiring circumstellar interaction. For CSM 
interaction models to be the sole explanation, they must 
have CSM mass comparable to the ejected mass (at a large 
radius). In addition they require the CSM to have a surpris¬ 
ingly narrow range of densities across all the objects. 

• The shape of the light curve around peak is intrinsically 
very similar to that of normal SNe Ibc, except SLSNe are 


11 CONCLUSIONS 

We have investigated the physical properties of the largest 
sample of super-luminous supernovae constructed to date, 
and have found that they appear to be closely linked with 
normal SNe Ic. All of the properties of the class can be 
explained by taking a normal SN Ic, increasing the ejected 
mass, and re-shocking the ejecta using a powerful central en¬ 
ergy input, such as the emission from magnetar spin-down. 
The key observables are a peak luminosity boosted by ~ 3 
magnitudes, a broader light curve from the large diffusion 
mass, and a flat velocity evolution caused by the hydro- 
dynamical impact of the additional energy source. At late 
times, the light curve shape deviates from that of a normal 
SN Ic, which follows the radioactive decay of 56 Co after a few 
diffusion times from peak. In super-luminous objects, the 
56 Co decay is masked by some other dominant heat source 
(e.g. magnetar). 

Other central engine models have been proposed, for 



























































26 Nicholl et al. 


example fallback accretion onto a newly-formed black hole 
( Dexter k, Kasen|2013 ), but it is magnetar models that have 
been explored most in the literature. The fundamental lower 
limit on a neutron star spin period is ~ 1 ms, and Galac¬ 
tic magnetars have magnetic fields < 10 15 G. None of the 
SLSNe yet discovered—neither the slowly declining objects, 
such as PTF12dam ( Nicholl et al.||2013 ), nor the brightest 
SLSNe, like CSS121015 ( Benetti et al.|2014 )— have emitted 
an integrated luminosity higher than that expected for mag- 
netars with spin periods in the range 1-10 ms and magnetic 
fields of a few times 10 14 G. It would be a problem for the 
magnetar model of SLSNe if an object was discovered that 
could only be fit with a sub-ms spin period or magnetic field 
> 10 15 G, but otherwise looked like a typical SLSN Ic. This 
could be one way to discriminate between competing models 
with different central engines. 

Another important test will be detailed calculations of 
synthetic spectra. So far, only Dessart et al. ( 2012 ) have pre¬ 
sented calculations of model SLSN spectra based on magne¬ 
tar radiation. Howell et al. (2013) presented parameterised 
models, putting a bright central energy source inside carbon- 
and oxygen-rich ejecta, such as might be expected from a 
stripped-envelope SN. In both cases, good matches were 
found to observational data. The main competing theory, 
that SLSNe are powered by hydrogen-poor circumstellar in¬ 
teraction, will also have to pass this test, but calculating 
such a spectrum is at the limit of current modelling capa¬ 
bilities. Another issue is that the interaction model has so 
many tuneable parameters that a wide range of models will 
need to be produced to compare with observations. However, 
light curve fitting shows that spectral modellers should focus 
on the regime with several solar masses of ejecta and CSM, 
with pcsM ~ 10 -12 gcm -3 . (We also note that interaction 
models provide a good description of the photometric and 
spectroscopic evolution of SLSNe II, which in some cases 
share characteristics with SLSNe Ic; Benetti et al. 2014). 
Based on our analysis here, we prefer magnetar-like models 
for SLSNe Ic, at least until such time as synthetic spectra 
exist for the CSM model. 

On the observational side, future work should focus on 
finding more low redshift SLSNe, in order to improve on the 
statistics presented here. This should reveal more clearly 
whether, for example, all SLSNe Ic are drawn from a con¬ 
tinuous population, if there is an excess of very massive ob¬ 
jects, etc. Understanding the properties of these remarkable 
events is essential, particularly in preparation for the com¬ 
ing era of the James Webb Space Telescope and 30m-class 
ground-based facilities, which will allow us to find SLSNe at 
redshifts up to z ~ 10. By constraining their physics at low 
redshift, we may be able to use SLSNe to probe early cos¬ 
mological expansion, high-redshift dwarf galaxies, and the 
first generation of star formation in the Universe. 
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